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SUMMARY 


The different types of general instability already discussed in 
literature are critically reviewed. The writer’s earlier theoretic 
calculations of the critical load of the inward bulge type of general 
instability are revised in the light of the experimental data of the 
GALCIT investigation. A new simple formula is derived which 
is in reasonable agreement with experiment. 


TYPES OF GENERAL INSTABILITY 


ENERAL INSTABILITY is defined as the simultaneous 
buckling of the longitudinal and the transverse 
reinforcing elements of a monocoque structure. The 
definition is meaningless unless the monocoque has both 
longitudinal and transverse reinforcements, and only 
such monocoques are considered in the present paper. 
For this reason the word ‘‘monocoque”’ can be used in 
place of the more cumbersome ‘‘semimonocoque”’ or 
“reinforced monocoque”’ without incurring any danger 
of misunderstanding. 

There are several known types of instability which 
conform with the above definition, and undoubtedly 
additional ones will be described in the future. They 
vary considerably in the value of the critical loads and 
in the pattern of the distortions. To facilitate an under- 
standing of the problems involved it appears advanta- 
geous to adopt a different name for each of the various 
types of general instability. Such names are suggested 
below for the types already described in literature, and 
their characteristic features are briefly discussed. 


Flattening in Compression 


Wang'* showed that a monocoque cylinder loaded 
with uniformly distributed compressive forces may fail 
by flattening. ‘‘Flattening’’ means that the originally 
circular rings (transverse reinforcing elements) distort 
to ellipses and the originally straight stringers (longi- 
tudinal reinforcing elements) to sine curves of one or 
more half-wave lengths (see Fig. 1). The critical load 
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Fic. 1. Flattening type of general instability. 

was calculated by the Rayleigh-Ritz-Timoshenko 
method. The distortions of the rings and stringers 
were assumed to be inextensional and were represented 
by trigonometric series. The sheet-metal covering was 
assumed to be infinitely thin and the strain energy stored 
in it was neglected. 

Experimental evidence of this type of general in- 
stability is not available. It might be obtained with 
monocoques that are so proportioned that the buckling 
stress of their panels of sheet is much smaller than the 
buckling stress corresponding to general instability, 
provided that the reduction in the resistance to shear 
does not result in some form of torsional instability. 
If, however, the panels of sheet are still in the un- 
buckled state when general instability should occur 
according to Wang’s calculations, the strain energy 
due to shear in the sheet covering is not negligibly small, 
and it is likely to influence greatly both the critical 
load and the distortion pattern. 


Diamond Pattern in Compression 


The deflections are sinusoidal in both the longitudinal 
and circumferential directions and usually consist of 
several half-waves in either direction (see Fig. 2). The 
pattern is much the same as the one well known from 
the buckling of thin-walled nonreinforced circular 
cylinders. By definition the half-wave length in the 
longitudinal direction must be greater than the spacing 
of the rings if the buckling is to be classified as general 
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Fic. 2. Diamond pattern type of general instability. 


instability. The flattening discussed previously is 
just that special case of the diamond pattern when the 
number of half-waves is four in the circumferential di- 
rection. 

The buckling load of the diamond pattern type of 
general instability was calculated by Fliigge,* Dschou,? 
Taylor,'! and Ryder.’ The problem is also discussed 
by Timoshenko.'* In all these calculations the actual 
monocoque is replaced by an equivalent orthotropic 
shell—that is, by one that has different but uniformly 
distributed bending rigidities in the longitudinal and 
circumferential directions. The shearing stresses in the 
sheet covering are not disregarded. 

Experiments carried out by Nissen® with monocoque 
cylinders covered with corrugated sheet showed good 
agreement with Dschou’s theoretic formulas if experi- 
mentally determined values were used for the stiffness. 

Many attempts were made to use the theoretic 
formulas derived for the diamond-type buckling in 
compression for the calculation of the critical stresses 
of monocoque cylinders tested in bending. The basis 
for such a procedure is the argument that the critical 
stress in bending must be substantially the same as that 
in compression if the wave formation occurs only in the 
close proximity of the most highly stressed compressive 
fiber, since then the stress is materially constant in the 
region involved. Critical stresses predicted by the 
theoretic formulas of the diamond pattern general 
instability were checked at the Guggenheim Aeronau- 
tics Laboratory of the California Institute of Tech- 
nology (GALCIT).‘ It was found that the theoretic 
critical stress was invariably higher than the experi- 
mental stress, in some cases more than ten times as 


much. 

Two reasons may be put forth in explanation of the 
discrepancy. First, it is likely that in most of the speci- 
mens tested at GALCIT the wave length in the circum- 
ferential direction was so large that the variation of the 
direct stress due to bending had a decisive influence 
upon both the distortion pattern and the critical stress. 
If this is true, the instability cannot be classified as of 
the “diamond pattern in compression” type and the 
critical stress has to be calculated as suggested in the 
last section of this paper. Second, in-all the test 
specimens the wave formation was restricted to areas 
where the sheet covering buckled long before general 
instability occurred. Thus the resistance of the sheet 





to shear was materially reduced. It is conceivable that 
a modification of Dschou’s or Taylor’s theoretic for- 
mulas to take care of this effect would result in better 
agreement between theory and experiment. 

One might be inclined to expect better agreement 
also with bent monocoque cylinders in which the sheet 
covering does not buckle before general instability oc- 
curs. This can be the case in cylinders of large radius 
reinforced with closely spaced weak stringers and rings. 
It is quite likely, however, that the “‘linearized”’ theory 
of elasticity does not suffice for the calculation of the 
critical stresses of such structures, since it is known from 
the investigations of von Karman, Tsien, and Dunn’ 
that a nonlinear theory is necessary to explain the dis- 
crepancy between theory and experiment in the case 
of the buckling of thin-walled nonreinforced cylin- 


ders. 


Flattening in Bending 

Brazier! showed that an isotropic thin-walled cir- 
cular cylinder flattens to an elliptic shape when a pure 
bending moment is applied to it. The flattening in- 
creases gradually when the applied moment is increased, 
and the cylinder collapses when the moment reaches its 
critical value. 

Brazier’s calculations were extended by Heck® to 
describe the behavior of orthotropic cylinders of elliptic 
cross section and by the writer (Part 2)° to describe the 
behavior of circular cylindrical structures consisting of 
stringers, rings, and a very thin sheet covering. Ac- 
cording to all these investigations, buckling occurs after 
the diameter of the cylinder foreshortens gradually to 
seven-ninths of its original length. This circumstance 
forcibly rules out the “flattening in bending” type of 
general instability as a form of buckling of actual air- 
craft structures. No designer or licensing authority 
is likely to consider an aircraft structure airworthy if its 
elastic distortions during actual flight attain such mag- 
nitude. On the other hand, in large airplanes it might 
be well to calculate, by the theoretic formulas of Part 2 
of reference 6, the amount of flattening, in order to 
check whether the structure has adequate stiffness, 
and the stresses caused by flattening, in order to ascer- 
tain that the structure has adequate strength. The ex- 
periment reported by the writer® showed reasonable 
agreement with calculations. 


Inward Bulge in Bending 

Because of the lack of a better expression this name 
has been chosen to indicate that the characteristic 
feature of this type of buckling is ati inward bulge that 
develops symmetrically to the most highly stressed 
fiber in compression (see Fig. 3). In some cases second- 
ary bulges also appear and the deflected shape slightly 
resembles the diamond pattern in compression already 
described. However, deflections are greatest near the 
most highly stressed fiber in compression, they decrease 
in magnitude toward the neutral axis, and the tension 





“<2 








® 
i 


INSTABILITY OF MONOCOQUE CYLINDERS 107 





Fic. 3. Inward bulge type of general instability. 


side of the cylinder remains smooth. This type of 
buckling was theoretically derived and in one single 
experiment was verified by the writer in 1937.° In the 
calculations the monocoque cylinder of circular cross 
section was assumed to consist of stringers, rings, and a 
very thin sheet covering. The buckling load was de- 
termined by the Rayleigh-Ritz-Timoshenko method, 
and the strain energy stored in the sheet covering was 
neglected. 

Experimental results of an extensive series of inves- 
tigations are available for checking the validity of the 
theory of the inward bulge type of buckling. These 
experiments were carried out at GALCIT as part of a 
research program undertaken to clarify the problems of 
general instability. Results were submitted to the 
Civil Aeronautics Authority in five reports between 
1938 and 1940.4 A summary of the investigations may 
also be found in Sechler and Dunn’s textbook.” It 
appears that at buckling most of the specimens showed 
the inward bulge pattern. Nevertheless, a comparison 
at GALCIT of the test results with the writer’s theo- 
retic predictions showed marked disagreement. In 
the following sections of the present paper the assump- 
tion that caused the disagreement is discarded and the 
theory is amended accordingly. The revised theory 
yields simple formulas that show reasonable agreement 
with experiment. 


CALCULATION OF THE CRITICAL LOAD OF THE INWARD 
BULGE TYPE OF GENERAL INSTABILITY 


The original calculations of the critical load of the 
inward bulge type of general instability were published 
by the writer in Part 3 of reference 6, where this type 
of instability was termed “local buckling of a mono- 
coque structure.”” The calculations deal with a mono- 
coque cylinder of circular cross section consisting of 
evenly spaced rings, evenly spaced stringers, and a 
very thin sheet covering. The cylinder is under the 
action of a constant bending moment applied at the 
end sections by means of loads distributed in conform- 
ity with the Bernoulli-Navier assumptions. It is 
further assumed that the structure distorts in accord- 
ance with the Bernoulli-Navier theory of bending until 
buckling occurs. 

The stability of this state of equilibrium in bending 
is investigated by the Rayleigh-Ritz-Timoshenko 


method by determining the change in the total potential 
energy of the system during a virtual displacement. 
The displacement pattern is arbitrarily assumed to be 
that shown in Figs. 3, 4, and 5. It is given analytically 























Fic. 5. Deflected shape of stringer according to assumptions. 


by the expressions for the radial and tangential dis- 
placements w, and w,, respectively, of any point on the 
surface of the cylinder: 
w, = —(a/2)[1 — cos (2rx/L)|[cos (ng) + l 
cos (2n¢)]} (1) 
(a/2)[1 — cos (2rx/L)|[(1/n) X 
sin (ng) + (1/2n) sin - 


Ws 


provided that 
gS a= n/n (la) 


These equations satisfy a number of boundary condi- 
tions and represent inextensional deformations of the 
rings and stringers. They contain two arbitrary para- 
meters—namely, L, which is the wave length in the 
longitudinal direction, and the number 2 = r/o, 
which characterizes the wave length in the circumferen- 
tial direction. Next, the strain energy due to the dis- 
tortions according to Eqs. (1) is calculated, as well as 
the work done by the applied loads. The expressions 
obtained do not have terms with a; this circumstance 
indicates that the original undisturbed shape is a state 
of equilibrium (principle of virtual displacements). 
Whether or not this equilibrium is stable depends in 
the present case upon the terms containing a®. The 
critical value of the bending moment is reached—that 
is, the change from stability to instability occurs— 
when the strain energy stored in the structure is equal 
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to the work done during the virtual displacement. 
From this condition the critical load is calculated: 


Pea = (xEI,,/L*)[4eMs + (1/2*)(1/y)Mi]/(Mi + 
Ri) (2) 


where 
x = (r/L)*(r/d) (Ese sn/NE,T,) (3) 
and 
Deir = Isr, + (5/8)(1/n?)I ser, (4) 
Furthermore, 


M, = sin (x/n){1 — [1/(1 — m*)] + [1/2(1 —- 
4n?)] — [1/(1 — 9n*)] + [1/2(1 — 16n?)]} 
Mz = 17n*? — 10n + (2/n) 
M; (1/n) (5) 
Ri = (1/2n?) sin (x/n){ (5/4) — [1/(1 — n?)] — 
[1/(1 — 4n?)] + [1/(1 — 9m*)] — [1/401 — 
16n?)] } 


Il 


In these equations 


P., = the critical load of the most highly stressed 
stringer in compression 


I, = the moment of inertia of a stringer for radial 
bending 
Is, = the moment of inertia of a stringer for tan- 
gential bending 
I, = the moment of inertia of a ring (for radial 
bending) 
Ey, = the modulus of the material of the stringer 
E, = the modulus of the material of a ring 
L = the wave length in the longitudinal direction 
ry = the radius of the cylinder 
d = the stringer spacing 
N = a number dependent on the number m of 


rings that distort within one wave length 
during buckling 


Values of N are given in a table® for m = 1, 2, 3, 4, 
and 5. 
With L, denoting the ring spacing one has 


L= (m+ 1)i (6) 


The magnitude of the critical load given in Eq. (2) 
depends upon the two parameters 7 and m = (L/L) — 
1. At buckling the parameters must have values that 
make P., a minimum. The critical load® is minimized 
in the following way: The differential coefficient of 
P., with respect to m is calculated and equated to zero, 
and between the equation so obtained and Eq. (2) y 
is eliminated. The critical load may then be given in 
the form 


Pa = WE gerd gtr/ (ALi)? (7) 


where } is the equivalent length coefficient. The value 
of \ depends upon the number m of rings involved in 
buckling and upon a structural coefficient A given by 
the equation 





A = (r°/Li5d) (Ese s:/E,,) (8) 
The minimum value of P,, with respect tom = (L/L,) 
— lis then determined graphically with the aid of Fig. 
19° in which the equivalent length coefficient \ is 
plotted against the structural coefficient A for different 
values of m. From this chart the buckling load of the 
most highly stressed stringer of any actual monocoque 
cylinder can be easily calculated. One has only to de- 
termine the structural coefficient A of the particular 
monocoque structure and take from the chart the cor- 
responding highest value of \ compatible with the 
number m of the rings which can be involved in buck- 
ling. 

A printing error in reference 6 should be pointed out 
here. The expression in the first parenthesis on the 
right-hand side of the last of Eqs. (5) of the present 
paper (Eq. 78-a°) was printed erroneously as ‘‘!/2n?”’ 
instead of the correct “‘1/2n*.”” Moreover, Eq. (87) of 
reference 6 was misinterpreted at GALCIT. Its 
right-hand side was originally printed as [1/(m + 1)*- 
N,,|A. It was read as {[1/(m + 1)5]N,,}A instead of 
the correct {1/[(m + 1)* N,,]}A. 


EFFECT OF THE SHEET COVERING 


Reference 4 contains several diagrams and photo- 
graphs showing the distortions after failure of the speci- 
mens tested at GALCIT. Some of the specimens dis- 
played in the figures are of the regular sheet-covered 
monocoque type; others have wire bracings across the 
panels in place of the sheet covering. The distortions 
of the latter show good agreement with the assumptions 
stated in the preceding section. The former group in- 
cludes some test pieces that failed definitely according 
to the inward bulge pattern; others, the buckling of 
which might possibly be classified as being of the dia- 
mond pattern type; and still others, in which the wave 
length in the longitudinal direction was equal to the 
ring spacing so that the buckling was clearly not of the 
general instability type. 

Although the distortion pattern is not necessarily 
the same at incipient buckling and after failure has oc- 
curred, one might expect reasonably good agreement 
between theory and experiment in the value of the 
buckling load if the distortions measured after failure 
were similar to those assumed in the calculations. 
Nevertheless, at GALCIT the buckling loads obtained 
by experiment were from one and two-tenths to nine 
times greater than those calculated by the theory. 
The reason for this discrepancy may be found in the 
effect of the sheet covering and of the wire bracing in 
the two kinds of specimens 

The buckled shape given by Eqs. (1) corresponds to 
inextensional (purely bending) deformations of the 
rings and stringers. The same distortion pattern, how- 
ever, involves shearing deformations of the covering 
which were disregarded in the theory when the sheet 
was assumed to be infinitely thin. If the covering is 
in the nonbuckled state, the strain energy due to the 
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shearing strain is many times greater than that due to 
the bending of the stringers and rings. On the other 
hand, the resistance to shearing deformations is ex- 
tremely small in the panels that had buckled before 
the occurrence of general instability. Thus it is reas- 
onable to assume that the inward bulge is restricted to 
portions of the cylinder where the panels of sheet cover- 
ing are in the buckled condition. This statement, how- 
ever, means that the wave length in the circumferential 
direction is determined by the condition of the panels 
of sheet and, consequently, must not be calculated from 
the minimum requirement stated in the preceding sec- 
tion. In this connection it is worth noting that the 
envelope of the curves of Fig. 19° corresponds to the 
deformations extending over the entire circumference 
(mn = 1) and that the envelope contains the values of 
to be used if the monocoque cylinder is reasonably 
long. 

The determination of the length of the arc of cir- 
cumference involved in buckling is not yet possible 
without recourse to experiment. Among other things 
the length is influenced by the critical stress of the 
panels, by the maximum stress in the cylinder, and by 
the shift of the neutral axis of the cylinder due to the 
ineffectiveness of the buckled sheet in carrying com- 
pression. It is reasonable to assume, however, that 
¢o = 7/n is equal to, or smaller than, 90°. 

Similar considerations may be offered in the case of 
the wire-braced specimens. On the tension side of the 
monocoque the crossed wires are under tension and 
prevent changes of the right angles subtended by 
stringers and rings. On the compression side, however, 
the wires slacken and permit small changes in the angles. 
The length of the arc of circle involved in the inward 
bulge type of buckling depends, therefore, upon the 
amount of distortion necessary for general instability 
as well as upon the initial tension in the wires. A purely 
theoretic calculation of these effects is not yet pos- 
sible. 

On the other hand, if the value of ” is known, the 
formulas developed in the preceding section may be 
simplified considerably. 


REVISED THEORY 


If 2 is a known constant, the buckling load P,, of 
the most highly stressed stringer is a function of the 
longitudinal wave length LZ only. In such a case every 
symbol in Eq. (2) is a constant except ZL and y. These 
two can be expressed in terms of (m + 1) if one makes 
use of Eqs. (3) and (6) and calculates the value of N in 
Eq. (3). N is the number by which the strain energy 
stored in a ring at the middle of the longitudinal wave 
length must be multiplied in order to obtain the total 
strain energy stored in the m rings within one wave 
length. The strain energy in any one ring is propor- 
tional to the square of the maximum deflection. The 
maximum deflection 8 of the ring in the middle of the 
wave is —2a as may be seen from the first of Eqs. (1). 
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The maximum deflection of the i‘ ring is 


B, = —all — cos [2ri/(m + 1)]} (9) 
Consequently, 


i=m 


N = (1/4) ¥ {1 — cos [2ni/(m+1)]}* (10) 


It is shown in the Appendix that, by carrying out the 
summation indicated on the right-hand side of Eq. 
(10), 


N = (3/8) (m + 1) ifm = er 
By definition (11) 
N=lifm=1 j 


With the aid éf Eqs. (3), (6), and (11) Eq. (2) can be 
transformed: 





P = TE seed str 1 { 4rM3 
7 Li (M, + Ri) \(m + 1)? 
2 
3M: tnt (12) 
8 7 A 


In this equation £,,, is the modulus for the material of 
the stringer, J,,, is given in Eq. (4), Z; is the distance 
between adjacent rings, A is the nondimensional struc- 
tural coefficient given in Eq. (8), and the values of M,, 
M2, M3, and R; can be computed from Eqs. (5). The 
critical load can be calculated, therefore, if the number 
m of the rings involved is known. In along monocoque 
cylinder m will have the value that makes P,, a mini- 
mum. Assuming arbitrarily that N is a continuous 
function of m and is represented by the first of Eqs. 
(11), the minimum value of P,, can be determined by 
differentiating Eq. (12) with respect to m and equating 
to zero the differential coefficient. The equation so ob- 
tained can be solved for m: 


m = (32/3)!/4 x(Ms/Mz)”4Av4 — 1 (13) 


This is the number of rings that distort within one wave 
length when the buckling load is a minimum. Eq. (13) 
can be written in the form 


m= wAY4 ~ om 1 


In Fig. 6 computed values of m are plotted against 
(1/n). The curve so obtained may be represented with 
an accuracy sufficient for practical purposes by the 
formula 


v = (2.84/n) 
Thus 
m = (2.84/n)AY¥4 — 1 (14) 
Eqs. (13) and (14) are approximate, since by definition 
m must be an integer. Nevertheless, substitution of 
the value of m from Eq. (13) into Eq. (12) should give 


the critical load P., sufficiently accurately if the number 
m of rings required for a minimum of P,, is available. 
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Fic. 6. Number m of rings in a wave when P,, is a minimum. 















































; Poe = wE gel ete/ (ALi)? (15) 
where the equivalent length coefficient \ has the value 

\ = RAYS (16) 
and the coefficient k can be determined from the equa- 
tion 


k? = (x/+/6)[(Mi + Ri)/-V/ MM] (17) 


In Fig. 7 computed values of k are plotted against 1/n. 
It can be seen from Fig. 7 that a good approximate 
formula for k is 


k=1/n (18) 


With the aid of Eqs. (16) and (18) Eq. (i5) may be 
rewritten in the form 


Pa = eB se] ser/LY2AY* 
Upon substitution of A from Eq. (8), Eq. (19) becomes 


2 d rV/ Esl trl, (20) 
Ly r? 


(19) 


Py =f 





According to Eq. (20) the maximum load that the most 
highly stressed stringer can carry is equal to the buck- 
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Fic. 7. Value of the equivalent length coefficient \ when Pe is 
a minimum, 
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ling load of a column, the length of which is equal to the 
radius of the cylinder and the bending rigidity of which 
is equal to the geometric mean of the bending rigidities 
of stringer and ring, provided that an end fixity factor 
C = n*/d/L, is used. In the expression for the end 
fixity coefficient, (d/L,) is the ratio of the stringer 
spacing to the ring spacing and m = w/g. This formula 
may be given in several alternate forms if desired. One 
of them is 


€or = n>? (pstrpr/T?) (E,t,/E sitet) v3 (21) 


where e¢,, is the critical strain in the most highly stressed 
fiber of the monocoque cylinder and the other symbols 
not previously defined are given by the formulas 


= Ts:r/Aser pr” 7” I,/A, t 22) 
= Agi,/d = A,/Ii } fis 
A,,, and A, are the cross-sectional areas of a stringer and 
a ring, respectively, including the effective width of 
sheet. 


2 
Prey 


tetr 


Los,, A 
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Fic. 8. Theoretic and experimental values of \ for wire-braced 
specimens (m = 3), 


COMPARISON WITH EXPERIMENT 


The formulas for the buckling load developed in the 
preceding section are easy to use if the wave length in 
the circumferential direction and, consequently, the 
value of m are known. For the time being this value 
has to be determined by experiment. An analysis of 
the 22 wire-braced test specimens of the GALCIT 
investigations yielded m = 3 as the best value. With 
this, k = 1/3 from Eq. (18) and Eq. (16) be- 
comes 


\ = (1/3)A”4 (23) 


A graphic representation of this formula on logarithmic 
paper is a straight line with a slope equal to one- 
fourth. It is shown in Fig. 8. The figure also con- 
tains values computed from the GALCIT experiments. 
The abscissa of each experimental point was calculated 
from Eq. (8); the ordinate, from Eq. (15) in which the 











neglected. 














experimentally determined value was substituted for 
P., and the equation was then solved for X. 

The experimental points are in fair agreement with 
the theoretic straight line. The maximum deviation 
to the safe side is 20 per cent; the average deviation, 
10 per cent. The maximum deviation to the unsafe 
side is 26 per cent; the average deviation, 12.5 per cent. 

In Fig. 9 theory and experiment are compared on 
the basis of the 35 sheet-covered GALCIT test speci- 
mens that failed in general instability. The value of 
n is again assumed to be 3. With sheet-covered speci- 
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Fic. 9. Theoretic and experimental values of \ for sheet-covered 
specimens (m = 3). 


mens the radial and tangential moments of inertia must 
be calculated with proper consideration of the effective 
width of sheet, and the difference between the location 
of the centroid of the section plus effective width and 
that of the centroid of the section alone must not be 
In the case of the stringers Marguerre’s 
cube root formula was used for the effective width of 
the flat sheet, and the influence of the curvature was 
taken into account according to Ebner’s suggestion. 
The maximum stress measured in the tests was taken 
for the edge stress. The pertinent formulas are dis- 
cussed in Part 1 of reference 6. They are given below 
in a slightly modified form and with the notation of the 
present paper. 


2w’ = 2w + (€curvea/€maz-)(d — 2w) (24a) 


2w = [egat/(€maz- — unveadl “YS (24b) 
€aar = 3.62(t/d)? (24c) 
€curved = 0.3(t/r) (24d) 


In the above formulas 2w’ is the effective width of the 
curved sheet; €mgr. the strain measured in the most 
highly stressed stringer on the compression side of the 
test specimen; d, the stringer spacing; ¢, the thickness 
of the sheet metal covering; and 7, the radius of the 
cylinder. In the case of the rings the effective width of 
sheet is assumed to be equal to the width of the ring 
section for reasons discussed in paragraph 22.6 Except 
for the above details, the procedure in calculating the 
experimental points is the same in Figs. 8 and 9. 
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The experimental points of Fig. 9 show the general 
trend of the theoretic straight line, but the scattering 
is considerable. The maximum deviation to the unsafe 
side is 19 per cent; to the safe side, 44 per cent. The 
average deviation to the unsafe side is 7 per cent; to the 
safe side, 22 per cent. According to these figures some 
of the specimens failed under much higher loads than 
those predicted by theory on the assumption of m = 3, 
since in the equation of the critical load the value of \ 
appears to the second power. 

In order to get an insight into the reasons for the dis- 
crepancy, the values of ” that make theory and experi- 
ment agree were calculated for each specimen. Since 
enters into the calculation of J,,, (see Eq. 4), it was nec- 
essary to proceed in the following way: 

Using subscripts 3 and p to indicate that m is as- 
sumed to have the value 3 and p, respectively, Eq. (3) 
may be rewritten in the form: 


Po = WE str stry/(AsL1)? 
Substituting for P,, the value P,,, obtained in experi- 
ment, 
As = [eB seLstry/PexpLr] V2 (25) 
This is the experimental value of \ provided that m = 3. 
If n = p, 
Ny = [4*BserDetry/ Peep Lr] ¥* 
Consequently, 
Ny = AslLetrp/Leen] ¥? 
Moreover, from Eq. (8) 
Ay = As(Lstrp/Istrs) 
According to Eqs. (16) and (18), theory requires that 
dX = (1/n)Av4 
Thus 
A, = p*h,* 
Substituting the above determined values, 
As = p*rs* [Lserp/T str] 


Substituting for J;,,, its value from Eq. (4) and trans- 
forming the expression so obtained, 


DT str + 27(5/8)L sen ‘ans (Asl str,/As*) = 0 (26) 


This is a quadratic in p*. I,,,, and I;,,,, are properties 
of the stringer section, 7; and A; can be calculated from 
Egs. (4) and (8) with » = 3, and A; is the “‘experimen- 
tal’”’ coefficient calculated from Eq. (25). Solution of 
Eq. (26) for p gives the value of m that makes Eq. (16) 
hold. 

For every sheet-covered specimen the value of m was 
determined according to the above outlined procedure. 
It was found that m varied between 2.12 and 5.59. A 
qualitative explanation for the differences in the values 
of » is not difficult to find. The inward bulge can de- 
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velop only if the distortions of the stringers and rings 
are not restrained materially by the panels of sheet. 
However, the resistance of the sheet to distortions de- 
pends upon two circumstances: first, whether it is in 
the buckled state and, second, how much movability 
results from its buckling. 

The length of the arc of the circumference along which 
the panels had buckled depends upon the ratio of the 
stress in the most highly stressed fiber of the cylinder 
to the buckling stress of a panel of sheet. Since in most 
cases the curvature effect is decisive, the buckling stress 
of the panel may be assumed to be proportional to 
Et/r, where E is the modulus of the material, ¢ is the 
thickness of the sheet, and r is the radius of the cylinder. 
Thus the ratio of the greatest stress to the buckling 
stress is proportional tO €mar.7/t, Where €mar is the 
greatest strain measured in the experiment. This 
quantity is one of the parameters that may be expected 
to have influence upon the value of n. 

The amount of movability resulting from buckling 
depends upon the ratio of the difference between the 
length of the arc of the panel and the length of its chord 
to the length of the chord. In a first approximation 
this ratio is equal to r/d, as can be proved without dif- 
ficulty. Consequently, r/d is the second parameter in 
question. 

In Fig. 10 calculated values of are plotted against 
the two parameters (€mgz.7/t) and (r/d). Curves of 
constant ” are drawn in a way to give the best agree- 
ment with experiment. The figure shows that n is 
large if €mar.7/t is small—that is, if the sheet had 
buckled along only a small part of the circumference— 
and if r/d is large—that is, if buckling caused only a 
small amount of movability. Moreover the effect of 
Emax.‘/t upon m is small if r/d is small—that is, if the 
movability of the buckled panels is great—and_ the 
effect of r/d upon n is small if énaz.7/t is large—that is, 
if the buckling stress is well exceeded in many panels of 
sheet. 
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Experimental values of m for sheet-covered specimens. 


It may be seen from Fig. 10 that the differences are 
small between the values of m according to the curves 
and according to the experimental points. The maxi- 
mum difference is 14 per cent and the average differ- 
ence is 4.4 per cent, if one does not include the four 
points for which m is equal to 5.03, 2.45, 2.37, and 2.12. 
With the first two, failures started by panel instability; 
with the last two, a consistent application of the formu- 
las resulted in an exaggerated increase in the value of the 
moment of inertia due to the effective width of sheet. 
In the case of the last two points the differences are: 
18 per cent for point 2.37 and 29 per cent for point 2.12. 

In all the above calculations it was assumed that the 
test cylinders were long enough to develop the distor- 
tion pattern corresponding to the minimum value of the 
critical load. The validity of this assumption was 
checked with the aid of Eq. (14). It was found that the 
number m of rings required for a minimum of the criti- 
cal load was available in all but the specimens plotted 
at the lower left corner of Fig. 10. A recalculation of 
the critical load of these specimens with the aid of Eq. 
(12) was not carried out, since the accuracy of the com- 
puted critical loads of these specimens is doubtful for 
the reason discussed under paragraph 5 in the following 
section. 


VALIDITY OF THE CALCULATIONS 


The validity of a theoretic calculation depends upon 
the assumptions made. In a discussion of the validity 
of the present calculations the following items must be 
considered : 

(1) Assumption of inextensional deformations. 

(2) Assumptions concerning the distortion pattern. 

(3) Assumptions concerning the circumferential 
wave length. 

(4) Method of calculating the buckling stress of a 
curved panel. 

(5) Method of calculating the effective width. 

(6) Assumption of the same moment of inertia for 
each stringer. 

(7) Assumption of the validity of Hooke’s law. 

(8) A disregard of the shift of the neutral axis. 


A discussion of these items follows. 

(1) It is believed that the assumption of inexten- 
sional deformations of stringers and rings is satisfac- 
tory if the panels of sheet covering buckle well before 
general instability occurs. In such a case the strain 
energy stored in the sheet due to the inward bulge type 
of instability is negligibly small, and the strain energy 
stored in the stringers and rings is a minimum for in- 
extensional distortions. If, however, general instability 
occurs before the buckling of the panels, the strain 
energy stored in the three elements—namely, stringer, 
ring, and sheet—is likely to be a minimum for a distor- 
tion pattern in which the distortions of stringers and 
rings are not inextensional. The present theory does not 
include this type of general instability, the calculation 
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of which would probably require the development of 
a nonlinear theory of the orthotropic cylinder along the 
lines of reference 7. , 

(2) It would be possible to make the assumptions 
concerning the distorted shape of the rings conform 
better with experiment by adopting a new expression 
for the distortions containing several parameters. It is 
likely, however, that the deviations of the present as- 
sumptions from the true deflected form have a smaller 
influence upon the critical loads than do the other items 
discussed below. Thus a revision of the assumptions 
concerning deflections does not appear warranted be- 
fore the other more important questions are clarified. 

(3) Because of the inaccuracy of the assumptions 
concerning the distortions near the boundary of the 
bulge (¢ = go) the value of m has no real physical sig- 
nificance. The angle gq = 7/n approximately repre- 
sents the limit of the major distortions and thus roughly 
coincides with the angle beyond which the panels of 
sheet are in the buckled state or beyond which the wire 
bracing is slack. The actual value of m, however, must 
be determined from experiment. It was found to be 
approximately 3 for the wire-braced specimens and to 
vary considerably in the case of the sheet-covered speci- 
mens, as shown in Fig. 10. 

(4) The buckling stress of the curved panels of 
sheet is of importance in the calculation of the effective 
width, which in turn influences materially the value of 
the critical load in general instability. Unfortunately, 
the actual value of the numerical coefficient in Eq. (24c) 
is not necessarily 0.3, and it may even vary from speci- 
men to specimen. However, because of the uncertainty 
of the basis of the theory of buckling of thin curved 
sheets and in view of the not clearly defined edge condi- 
tions, a refinement of the calculations involving the use 
of different coefficients in different specimens does not 
seem to be warranted. 

(5) The calculation of the effective width of sheet 
with the aid of Eqs. (24a) and (24b) is not reliable, but 
no better method is known to the writer. General 
practice is followed in the use of the effective width so 
obtained in the calculation of the radial moment of 
inertia of the stringer plus effective width combinations. 
On the other hand, the adoption of the same effective 
width for the calculation of the tangential moment of 
inertia was entirely arbitrary. This circumstance is be- 
lieved to have led to exaggerated values of J; and, 
consequently, to too low values of m in the case of the 
specimens marked 2.37, 2.12, 2.68, 2.82, 2.86, and 2.63 
in Fig. 10. 

From a theoretic standpoint the reduced effective 
width should be used in the calculations instead of the 
effective width. This means that in Eq. (24b) the 
right-hand side should be multiplied by two-thirds. 
The factor was omitted, since the average effective 
width over all the stringers involved in general insta- 
bility is greater than the effective width of the most 
highly stressed stringer, as discussed under (6). 


INSTABILITY OF MONOCOQUE CYLINDERS 
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(6) In the derivations it was assumed that the 
moment of inertia J,,, was the same for every stringer 
involved in buckling. In reality this moment of inertia 
varies with the location of the stringer. Thus some 
average value should be used rather than the minimum 
corresponding to the most highly stressed stringer, but 
in what way this average should be calculated is un- 
known. The procedure of calculating J,,, for the most 
highly stressed stringer was adopted for two reasons: 
first, it gives a clearly defined rule for the computations 
and, second, the most highly stressed stringer contri- 
butes the greatest amount to the expressions in the 
energy calculations. 

(7) It was assumed that Hooke’s law was valid for 
all the elements involved. In some of the test pieces the 
maximum strain €m,,, was above the proportional limit. 
Since the stress-strain curve of the material was not 
known, it was impossible to determine accurately the 
effect of exceeding the proportional limit. It is not 
likely, however, that this circumstance would apprecia- 
bly influence the results. 

(8) The neutral plane of the bent cylinder was as- 
sumed to be the plane of symmetry of the unloaded 
cylinder. This is not true when part of the sheet cover- 
ing is ineffective due to buckling, but this circumstance 
is not likely to influence appreciably the results. 

The above discussed shortcomings of the theory are 
balanced by the use of a value m calculated from ex- 
periment. Thus it is reasonable to believe that the 
critical load in the inward bulge type of instability of 
monocoque cylinders having proportions similar to the 
GALCIT test specimens can be predicted with an ac- 
curacy sufficient for practical purposes if the procedure 
outlined in the next section is followed. For monocoque 
cylinders that differ materially from the GALCIT 
specimens no such definite statement can be made with- 
out further tests. In many practical cases, however, a 
reasonable amount of agreement may be expected, since 
the present procedure seems to work for such dissimilar 
structures as the wire-braced cylinders and the sheet- 
covered monocoques having rings of small cross section. 

If some of the above-listed assumptions were re- 
placed by different ones, the curves of Fig. 10 would 
have to be modified. For instance, a change in the cal- 
culation of the tangential moment of inertia would ma- 
terially alter the curves in the lower left corner but 
would influence the rest of the figure only slightly. It 
is hoped that future investigations will succeed in es- 
tablishing a more rational basis for the determination 
of the coefficient 1. 

From a comparison of the formula developed in the 
present paper with the empirical GALCIT formula of 
page 333, reference 10, the following conclusions may 
be drawn: 

The GALCIT formula gives excellent agreement in 
the case of the test specimens that have the smallest 
ring section. In these the ratio of the moment of inertia 
of the stringer to that of the ring is 24.3 if the effective 
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width of the sheet is disregarded. For the calculation 
of the instability of these test pieces the GALCIT 
formula is superior to the formula of the present paper. 

In the case of the five test specimens having the two 
heavier ring sections the ratio is 1.36 and 1.20. With 
these the agreement of the GALCIT formula is not 
quite so good. Critical loads calculated by it are just 
about as close to the experimental values as those cal- 
culated by the formula of the present paper. This is 
true if the curves of reference 4 are considered; in Fig. 
8-13 of reference 10 two of the points are replotted on 
the assumption of a reduced effective radius of gyra- 
tion. Theory does not appear to give a reason for such 
a reduction, since the rings are not likely to be under 
the action of high compressive stresses before general 
instability occurs. Nevertheless, the reduction would 
be perfectly acceptable if supported by sufficient experi- 
mental evidence. This, however, is not yet avail- 
able. 

In the case of the wire-braced specimens the GALCIT 
formula fails to give agreement with experiment. 
Sechler and Dunn do not suggest its use for this type of 
structure, and the experimental points are not included 
in Fig. 8.13 of reference 10. On the other hand, the 
formula of the present paper leads to reasonable agree- 
ment with experiment in the case of the wire-braced 
specimens. Consequently, the predicted critical load 
of a monocoque cylinder in which the moment of inertia 
of a ring is greater than the moment of inertia of a 
stringer should be regarded with greater caution if 
calculated by the GALCIT formula than if calculated 
by the procedure of the present paper. 


PROCEDURE SUGGESTED FOR THE CALCULATION OF THE 
CRITICAL STRESS OF THE INWARD BULGE TYPE OF 
INSTABILITY IN PURE BENDING 


(1) Assume the value of the critical strain é€mer 

(2) Calculate the effective width 2w’ of the sheet by 
Eqs. (24). 

(3) Determine the moments of inertia of the stringer 
plus effective width for bending in the radial direction 
(Ist) and in the tangential direction (J;,,,). Calcflate 
(I;:,) with the aid of Eq. (4). Assume = 3 if no more 
accurate value is available. 

(4) Calculate the moment of inertia J, of the ring, 
assuming that the effective width of sheet is equal to 
the width of the ring. 

(5) Calculate the critical load P,, of the most highly 
stressed stringer from Eq. (20), assuming » = 3 if no 
more accurate value is available. The maximum com- 
pressive stress is then P,,/A;,,, where A,,, is the cross- 
sectional area of the stringer plus the effective width 
of sheet and ém.; is the maximum stress divided 
by E. 

(6) If €mar, calculated according to (5) differs 
greatly from its value assumed under (1) and if the cor- 
responding value of m in Fig. 10 differs greatly from its 
value assumed under (3) and (5), repeat the procedure 
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with new assumptions for €,.;, and until satisfactory 
agreement is reached. 

Eq. (14) may be used for ascertaining that the num- 
ber m of rings that distort within one wave length, 
when P,, is a minimum, is not greater than the number 
of rings available. If it is greater, Eq. (12) should be 
used for calculating the critical load, and the number of 
rings available should be substituted for m. It should 
be noted, however, that the use of Eq. (20) instead of 
Eq. (12) causes an error on the safe side. 

If the moment of inertia of the ring is equal to, or 
greater than, the moment of inertia of the stringer and 
if the sheet covering is thin, the assumption of m = 3 is 
likely to give a reasonable approximation. With this 
assumption Eq. (20) becomes 


Py = sssaq/t VE seliyE rl, 


: (27) 
1 


Tete = Ietre + 0.06944I yr (28) 


ACKNOWLEDGMENT 


The writer is indebted to Eugene S. Davidson for 
the analysis of the experimental data and the compari- 
son with theory. This work was continued by Conrad 
Crown, who also made theoretic calculations, in par- 
ticular those contained in the Appendix. 


APPENDIX 


The summation indicated in Eq. (10) can be carried out in the 


following way: 
2xi |? 
m+1 


2at + 12) (14 4ni 

m-+1 _ m+1 
(3/8)m — ~ oe yt a 
ms og 3 ek 


i=1 t=1 
The above two sums can be calculated if the following transforma- 
tion is considered: 
F = cosy + cos 2y + cos 3y +... + cos (m — 1)y + cos ny = 
{1/(2 sin y)][2 sin ycos y + 2sinycos2y + 2sinycos3y+ ... 
+ 2 sin y cos (n — 1)y + 2 sin y cos ny] 


2 sin A cos B = sin (A + B) + sin (A — B) 
the following formula is obtained: 
F = [1/(2 sin y)][sin 2y + sin 3y — sin y + sin 4y — sin 2y + 
... + sin my — sin (n — 2)y + sin (wm + 1l)y — sin (m — 1)y] 
In the brackets all but three terms cancel out. Thus 


(Continued on page 130) 

















A Tabular Method of Propeller Blade 
Stress Analysis 





JOSEPH STUART, III* 


SUMMARY 


A tabular, trial-and-error method of evaluating the steady 
bending stresses in a rotating propeller blade is developed and 
examples of its application are given. The bending moment at a 
shank section of the blade is estimated and, by tabular integra- 
tion, the corresponding radial bending moment distribution is 
calculated. Additional trials with appropriately revised shank 
section bending moment estimates are made until the correct 
distribution is found, having the physically required moment 
value of zero at the blade tip. 

& Tilted blades and blades whose section centers of gravity do 
not in general fall on the blade axis are shown to be readily ana- 
lyzed. Because the final analysis can be quickly checked on a 
calculating machine and because the method is easily learned, 
versatile, and fast, its general adoption by the propeller industry 
is advocated. 


SYMBOLS 


R = propeller tip radius, in. 

radial distance to blade section from propeller axis of 
rotation, in. 

distance a point is forward of a reference plane per- 
pendicular to the propeller axis of rotation, in. 

= distance blade section c.g. is forward of the blade axis, 

in. 

= thrust per blade, Ibs. 

= thrust shear acting, Ibs. (= — 7" dT) 

centrifugal force acting through the blade section, Ibs. 

net bending moment about section c.g., in.Ibs. 

= moment of inertia of blade section about line passing 
through section c.g. and parallel to the chord, in.* 

blade material cross-sectional area, sq.in. 

modulus of elasticity of blade material, lbs. per sq.in. 

radial distance separating two adjacent moment or 
two adjacent slope stations, in. 
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INTRODUCTION 


7 IMPORTANT STEADY STRESSES in a rotating pro- 
peller blade are: (1) the direct tensile stress due to 
centrifugal force, and (2) the bending stress resulting 
from the net bending moment acting. 

The centrifugal stress (= C.F./A) is readily evalu- 
ated. The bending stress is readily evaluated when the 
net bending moment is known. 


DEVELOPMENT OF EQUATIONS 


Fig. 1 shows schematically a portion of a rotating 
propeller blade. As is customary in propeller blade 
Presented at the Power Plants and Propellers Session, Eleventh 
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Fic. 1. Propeller blade portion. 
stress analyses, the blade has been untwisted, set at 
zero angle, and loaded with a thrust loading only.' 
The blade length has been divided by primary, num- 
bered, moment stations into intervals X in. in length. 
The secondary, lettered, slope stations, also X in. apart, 
fall at the midpoints of the moment station intervals. 
Signs follow the conventions given in Niles and New- 
ell.? 

Slope, dz/dr, values always refer to the blade axis 
and never to the line passing through the blade section 
c.g.’s. The fundamental expressions for this stress 
analysis method are developed from elementary beam 
theory as follows: 


[oun dr + (dz/dr), = (dz/dr), (1) 
or, closely, 
X(M/ED); + (dz/dr), = (dz/dr), (2) 


Eq. (2) gives an adequately accurate value for the 
axis slope at b in terms of the axis slope at a and the 
bending moment at 1, the midpoint of the slope sta- 
tions’ interval. For future convenience, Eq. (2) is 
multiplied through by X and rearranged to give: 


[X2/(ED,:) Mi + X(dz/dr), = X(dz/dr), (3) 
The expression for the c.g. offsets of the running pro- 


peller is: 


[asian dr + (Z, — Z;) = (@ — a) (4) 


or, adequately closely, 


X (dz/dr), + (Z2 — Z;) = (% — %) (5) 
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In Fig. 1, Zz, and Z; are, of course, exaggerated. The 
expression for the moment is found to be: 


[rords+ [sar +a 


or, closely, 
(22 — a)(C.F.), + XS, + M, = 


The first term in Eq. (7) is the centrifugal restoring 
moment increment. 

As, EI, C.F., S, Z, and X values are known, Eqs. 
(3), (5), and (7) used in sequence give the slope at 3, 
(dz/dr),, and the bending moment at 2, M2, corre- 
sponding to the initial values (dz/dr), and M;. These 
equations permit the calculation of the moment dis- 
tribution of a blade from an assumed root station bend- 
ing moment value and known root slope, station by 
station out to the tip of the blade. 


M: (6) 


M, (7) 


TABULAR APPLICATION OF EQUATIONS 


The equations just developed lend themselves readily 
to tabular calculations. Fig. 2(A) shows the propeller 
blade selected for the examples, under static condi- 
tions. Its axis is then a straight line inclined forward 


at the tilt angle 6. The primary moment stations are 
designated by the solid transverse lines at 6-in. inter- 


vals, beginning with the 9-in. station and concluding at 
the 69-in. station, the tip station. The secondary 
slope stations are designated by the light, dashed trans- 
verse lines through the midpoints of the moment 
station intervals. 















































Fic. 2. Example propeller blade. 


Fig. 2(B) shows the same propeller blade rotating. 
Its axis is curved forward and its slope is now equal to 6 
only at the constrained shank section. The normal, 
positive thrust loading shown gives negative, thrust 
shear (S) values, and positive bending moments, tend- 
ing to place the camber side of the blade in compression. 
The centrifugal force, shear, J, and Z values of the ex- 
ample blade are plotted in Fig. 3. These values are 
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r, INCHES 
Fic. 3. Example blade basic data. 


assumed to have been previously obtained by known 
methods. 

Table 1 gives the final balanced trial for the example 
blade, for X = 6 in., E = 29,000,000 and 6 = 0.3° 
or 0.005233 rad. at the constrained, 6-in., initial slope 
station. In Table 1 values from Fig. 3 are multiplied 
or subtracted, as indicated, to give columns 3, 6, 8, and 
10 directly. 

The initial entry in Table 1 is the known slope entry 
underlined at the head of column 5: 


6 X 0.005233 = 0.03140. 


The other underlined entry is the 9-in., initial moment 
station, assumed root moment value of 28,032 in.Ibs. 
From these starting values and the given values in 
columns 3, 6, 8, and 10, Eqs. (3), (5), and (7) are ap- 
plied in sequence as follows: 


0.00355 + 0.03140 = 0.03495 


These initial entries in columns 4 and 5 are in ac- 
cordance with Eq. (38). 


0.03495 + (—0.169) = —0.13405 


These entries in columns 5, 6, and 7 are in accordance 
with Eq. (5). 
(—9,920) + (—7,680) + 28,032 = 


These entries in columns 9, 10, and 2 are in accord- 
ance with Eq. (7). 

These sequences are repeated station by station out 
to the tip moment station. In this particular example 
the tip moment has the physically required value of 
zero and the bending moments listed in column 2 de- 
fine the actual blade bending moment distribution. 

Table 2 gives an adequately balanced moment dis- 
tribution for the same blade analyzed in Table 1 save 


10,432 





PROPELLER BLADE STRESS ANALYSIS 





TABLE 1 





(1) (2) (3) (4) (5) 


36 36M dz 
— X 108 ar of 


M, In.Lbs. EI 
0.03140 


Station 


28,032 0.1267 0.00355 


0.03495 


10,432 0.1267 0.00132 
0.03627 


17,142 0.4320 0.00741 
0.04368 


11,601 0.9930 0.01152 

0.05526 
7,545 1.855 0.01400 

0.06290 


4,616 3.268 0.01509 
0.08429 


2,444 5.214 0.01274 
0.09703 

804 8.271 0.00665 
0.10368 

—0.00485 
0.09883 

—0.01801 


0.08082 


(6) 


Z2,—- Z1 


—0.13405 
0.21827 
0.03068 
0.05520 
0.07520 
0.08829 
0.09403 
0.09168 
0.08983 


0.29382 





that the blade section c.g.’s are now considered to lie 
along the blade axis (Z = 0) and the root tilt value has 
been taken to be zero. Comparison with Table 1 indi- 
cates that for blades without c.g. offsets two columns 
fewer are required. The facility with which a new 
blade tilt can be analyzed is worthy of note. 

For smaller blades or blades having large, erratic 
c.g. position shifts it may be desirable to take X equal to 
3 in. rather than equal to 6 in. The basic equations 
indicate that in Table 1, ‘‘9’’ would replace ‘‘36”’ in the 
headings of columns 3 and 4 and “3” would replace 
“6” in the headings of columns 5 and 10. The stations 
would read (7.5), 9, (10.5), 12, ete. 

These changes and their corresponding entry changes 
are the only complications entailed in changing X. 
However, to facilitate checking, it is recommended 
that X values of 6 or 3 in. be considered preferred 
values. 

To facilitate checking it is desirable that final trials 
be run on a calculating machine. The last significant 
figure should be carefully rounded off to permit identical 
figures to be calculated by the checker. 


ESTIMATION OF TRIAL RooT MOMENTS 


Usually only three trials are required to achieve a 
satisfactorily balanced moment distribution if the trial 
root moments are estimated as follows: 

(1) Take the first trial roof*moment estimate equal 
to 70 per cent of the moment the thrust would cause if 
it acted upon a stationary blade. 

(2) Take the second trial root moment estimate 


TaBLe 2 





(2) (4) (5) 


3 (# 
dr 


36 36M 
* x 10 
(g2 — 21) 


Station In.Lbs. 
(6) 0 

9 35,880 0.00455 
(12) 0.00455 

0.00362 
0.00817 

0.00926 
0.01743 

0.01502 
0.03425 

0.01840 
0.05085 

0.01944 
0.07029 

0.01648 
0.08677 

0.01101 
0.09778 

0.00434 
0.10212 


—7,680 


28,537 


537 —7,630 


21,444 


1,016 —7,330 


15,130 


1,648 —6,860 


9,918 


2,171 —6,140 


5,949 


2,411 —85,200 


3,160 


2,256 —4,085 


1,331 


1,760 —2,790 


301 
1,082 —1,375 
8 0.00022 


0.10234 348 — 346 





equal to the first estimate minus one-fourth of the tip 
moment found on the first trial. 

(3) Linearly interpolate, or extrapolate, results of 
the first and second trials, accurately, to estimate the 
root moment for the third trial. 

Fig. 4 shows a plot of tip moments against assumed 


root moments for the case considered in Table 1. The 
exact linearity displayed makes possible accurate inter- 
polation over great ranges. The origin of the one- 
fourth rule is evident in the indicated slope of this 
particular curve. 





118 





9" STATION MOMENT, IN.LBS. 


Fic. 4. Calculated tip moment versus assumed root moment. 

That only roughly one-fourth of an applied tip 
moment would appear at the blade root seems reason- 
able enough when one reflects that the application of 
the moment would cause deflections generating neu- 
tralizing restoring moments. 


SMALL UNBALANCED TIP MOMENTS 


It will be noted that Table 2 was not perfectly ad- 
justed but that an unbalanced tip moment of 10 in.Ibs. 
was allowed to remain. Fig. 5 shows the bending stress 
distributions calculated for another propeller blade 
with substantially larger unbalanced tip moments. 
From Fig. 5 it is concluded that small residual tip 
moments have no important effect upon stress values 
determined for the highly stressed portions of the pro- 
peller blade. It is recommended that moment distribu- 





Fic. 5. Effect of unbalanced tip moments on blade bending 
stresses. 
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tions for propeller blades, 10 ft. in diameter or larger, be 
considered adequately balanced when the residual tip 
moment is less than 50 in.Ibs. 

The above-noted insensitiveness of the blade bending 
stresses to small residual tip moments justifies the occa- 
sionally convenient practice of letting the final tip 
moment station fall as much as an inch beyond or an 
inch short of the actual tip of the blade. 

As a direct check of the method described, moment 
distributions were tabularly determined for a repre- 
sentative propeller blade under the following operating 
conditions: 

(1) Zero root tilt, take-off static thrust. 

(2) 0.5 degrees root tilt, take-off static thrust. 

(3) 0.5 degrees root tilt, windmilling at take-off 
r.p.m., with 10 per cent of maximum static thrust ap- 
plied as negative thrust. 

The tabularly determined moment distribution 
values were used as a basis for graphically determining 
deflections and corresponding centrifugal restoring and 
check net bending moment values. A comparison of 
the check and tabular distributions indicated that the 
greatest stress calculation error that would be encount- 
ered using the tabular method would be no greater than 
250 Ibs. per sq.in.—i.e., less than 1 per cent of the total 
steady stresses. 

As a check of the minimum number of moment sta- 
tions that it is safe to employ, the case given in Table 2 
was analyzed using 21-, 7-, 6-, 5-, and 4-moment sta- 
tions in addition to Table 2’s 11-moment station analy- 
sis. Agreement between the distributions was entirely 
satisfactory for the 7-, 11-, and 2l-moment station 
analyses. 

Corrected root bending moments checked, for all 
moment station numbers, to within 0.5 per cent. Tip 
moments checked less well, however, when fewer than 
7 stations were used. It is recommended that no fewer 
than 10-moment stations be used in final analyses. 


ADVANTAGES AND RECOMMENDATION 


In the two years that modifications of the method 
presented have been in use, three primary advantages 
have become evident: 

(1) The simple, repeated, tabular calculation se- 
quences are readily learned. 

(2) The numerically recorded tabular calculations 
are quickly checked on a calculating machine. 

(3) The method is several times faster than previous 
graphical methods. 

It is hoped that the method presented is sufficiently 
versatile and offers enough real advantages to warrant 
its early general adoption by the aircraft propeller 
industry. 
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ABSTRACT 


The analysis of bulkhead rings is at best a tedious process, 
and it is felt that any simplifications will be of help to the aircraft 
structural engineer. 

Of the several methods of analysis available, that outlined in 
reference 1 seems to be the simplest and most widely used. This 
paper is based on the general method of the N.A.C.A. report and 
presents, in particular, the application of the report to elliptical 
and circular frames. In addition, some general simplifications 
based on symmetry and load conditions are presented. 

No attempt has been made to break down the analysis into a 
number of specific solutions, but rather simplifications of the 
basic equations are presented which apply to any loading condi- 
tion. 

A set of curves is presented at the end of this article from which 
the length, moment of length, and moment of inertia of the length 
of an ellipse may be determined. 


INTRODUCTION 


— GENERAL EQUATIONS of reference 1 are modified 
in this paper to facilitate their application to sym- 
metric rings and rings of circular or elliptical shape. 
The modifications and simplifications are listed below. 
Steps I, II, III, and IV deal with specific properties of 
the sections and are progressive in that Step II em- 
bodies Step I, and Step III embodies Steps I and II. 
Steps V through VII depend on the type of loading. 


Step I—Simplifications Due to Symmetry. 

Step II—Cancellation of EJ Terms (Homogeneous 
Material—Constant Cross Section). 

Step III—Simplifications Due to Elliptical Shape. 

Step IV—Modifications Due to Splitting of Ring. 

Step V—General Loading (Unsymmetric). 

Step VI—Equal, Symmetrically Placed Loads. 

Step VII—Equal, Asymmetrically Placed Loads. 





Terminology 


Xm = bending moment at cut 

X, = axial load at cut 

X, = shear at cut 

M, = bending moment at any point in the deter- 
minate structure 

M, = bending moment at any point due to a unit 
bending moment at the cut 

M, = bending moment at any point due to a unit 
axial load at the cut 


Received April 21, 1942. 
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This figure illustrates the sign conventions used through- 
out this paper. 
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M, = bending moment at any point due to a unit Mr = Mog + MmgXm + MogXv + MopX, 
shear load at the cut lie 
ds = element of length waene, Oy eemaition: 
° . : Mm, = Mme = ] 
Subscripts Z and R refer to left- and right-hand sides 
° ° Moz = Mor = 7 ™ 
of the ring, respectively. see Fig. 1 
Moy, = & 
M, = Mo, + Mm,Xm + Mo,X> + Moy Xo Mop = =—% 


Step I 


The unknown loads at an assumed cut in the top of a symmetric ring are given below. The derivation of these 
expressions may be found in reference 1. 

Although the origin or point of action of these loads could have been assumed to be at any point in the plane 
of the bulkhead, past experience with bulkhead analysis has shown that from the point of view of work involved, 
in order to arrive finally at the stresses at any particular section of the bulkhead, the most logical origin is at the 
top or bottom of the bulkhead. It is to be remembered that the main part of the work in bulkhead analysis is 
not the determination of the unknown forces required to make the structure determinate but the determination 
of the effects of these redundant loads on the structure in connection with the known external loads and reactions. 
The assumption of an origin at the top of the ring leads to a simple convention of signs (as shown in Fig. 1). In 
addition, if the loading is complicated and a graphic method is resorted to for the determination of Mo, Vo, and 
Po, the calculations are simplified by the fact that each half of the frame can be treated as a separate cantilever. 


Of LS Ge) «+ J B)4)- 0/84) 4+ S)4] 
‘ 2S) (2s ar) - 2 Sz) 
(SNS (e+ Sen] 0S LS Ga), # + J GH), 


























Xp D) 

(ff a) - (fz) : 

X, = — [Sf (Mox/ED)z ds — Sf) (Mox/EID)p ds] /(2 f-x* ds/EI) (3) 
Step II 


If the ring is made of homogeneous material and is of constant cross section (or nearly so), the ‘EJ’ terms can- 
cel out and the above equations become: 








— (2f'y ds)[f (Moy), ds + J (Moy)e ds] — (2S y? ds)[S(Mo)z ds + S(Mo)e ds] (1a) 
4(f'ds)(f-y? ds) — (2f°y ds)? 
Wi te (2f-y ds)[f (Mo)z ds + {(Mo)eds] — (2f°ds) [Sf (Moy)z ds + S (Moy)e ds] (2a) 
4(f'ds)(Sy*ds) — (2f'y ds)? 
= —[f (Mor) ds — {(Mox)p ds] /2 fx? ds (3a) 
Srep ILI Evaluation of Integrals 
Eggs. (la), (2a), and (3a) are the general equations for Let 
a symmetric ring of any geometric shape. Since this Q, = complete elliptic integral of the first kind 
paper is concerned primarily with rings of elliptical or -/2 adden 
circular section, the integrals depending only on the = So" d6/V/) — k sin ?0 
geometric properties of the ring are completely deter- Q: = complete elliptic integral of the second kind 


mined below. 
The five basic integrals to evaluate are: (a) fo°/7ds, Se?V 1 — k* sin? 6 d0 


(b) So” y ds, (c)* So"? x ds, (d) So"? y* ds, and (ec) gx = Ge — 2) /p2 


3/2 
So"* x? ds. 
: , where a and 6 are the semimajor and semiminor axes 
of the ellipse, respectively, and @ is the basic parameter 


ll 


* Due to symmetry, this integral does not appear, but its 
evaluation proves useful for the solution of other problems. as shown in Fig. 2 
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(a) fo? ds = 2bQ. (b) SU? yds = 2670, ee ecancses~ 
(c) av x ds = ab/k(sin*k + k+/] — k?) 


's/2 


(d) f y? ds = 2b°Q, + 
0 


ons| (2 — ‘) , (C*) | 
- ( ae) @ + ap)? 
s/2 
(e) f x? ds = 2a7bQ. — 
0 


| (282 -— 1 1 — 
2a | ( 3k? ) a+ ( 3k a] 


Numerical values of these integrals are determined 
from graphs or tables at the end of this paper. It 
should be noted that the denominators and part of the 
numerators of Eqs. (la), (2a), and (3a) are now com- 
pletely determined. 











Strep IV 





(a) In many practical aircraft structures, it is neces- 
sary to make a cut in the ring. It is customary, when 
this is done, to run stringers from the cut to the ad- Fic. 3. This figure eat ope mage of X, for a split ring 
jacent bulkheads. , =F 

In the event that there are no stringers, the problem 
is no longer statically indeterminate because X,,, X,, and 
X, may be assumed to be zero at the free ends. 

Assuming now that stringers are present, let F = action Of Xm,Xp,X, 
resisting force offered by the stringer. It is assumed X es . x 
that the stringer has no torsional stiffness and that the 
stringer stiffness is negligible compared to the skin 
stiffness. 

It follows that under these circumstances X,, and 
X, will be zero at the cutout and that the stiffener de- 


Assumed poirt of 

















flection perpendicular to the plane of the skin at the Sy 

cutout is the only deflection that need be considered in | 

the equations of internal strain energy. 7 — 
@ 


Determination of Internal Energy of Stiffener 


The stringer is considered as a fixed-ended beam. 
Let: 








L = length of stiffener joining bulkheads adjacent 
to split ring 
¢ = distance from one bulkhead to split ring S/2 
d=L-e 
Then the deflection 6 = (1/3)(Fc%d*/EI,L'). 
U = internal energy of stiffener = (1/6) F*c*d*/EI,L* Fic. 4. Location of axes for split ring with fixed ends. 


dU/dF = (1/3) Fe*d*/EI,L* 


It should be noted that F = X,, since it has been as- 
sumed that this is the only force acting at the cut. 


r 373 Mox), ds — for x) p ds 2X, fea 
= - {# M,ds +**" =0 Ji " snaps sit 
OX, EI 3EI,L’ X.cidI X,cid*I 
; + = 0 
3L I, /t 3L°J, /Rr 


Since the ring is symmetric, 





“in 
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where J, = moment of inertia of the stringer and 
IJ = moment of inertia of section through side of 
ring. 


_ S (Mex), ds — f{ (Mox)p ds 
2S x? ds + (2/3)cdI/L*I, 





X, = 


The above equation assumes that the loads are dis- 
tributed in such a manner as to cause equal X, loads on 


either side. If the cut is small, this will normally be the 


case. 


Split Ring with Ends Rigidly Fixed 


(b) Bulkheads at wing spars are often discontinuous 
through the spar and rigidly attached to it. If the 
spar is assumed infinitely stiff compared to the bulk- 
head ring, no deflection can occur in the spar itself 
and no work can be done. If this is the case, the bulk- 
head becomes a split ring with unknown loads at the 
ends. 

A solution may be accomplished by making a hypo- 
thetic cut through the spar and integrating only over the 
periphery of the ring (see Fig. 4). It is necessary to 
integrate only over the ring itself, since it has been 
assumed that the spar is infinitely stiff compared with 
the ring. This fact may readily be detected from Eqs. 
(1), (2), and (3) where the stiffness is represented by the 
factor EJ, which in all the integrals appears in the 
denominator. It is obvious that if the EJ term over 
part of the path of integration is large compared to its 
magnitude over the rest of the path, then the contribu- 
tion of the first part to the total value of the integral is 
small, and no serious error is introduced by omitting 
the first part of the integration. 

Keeping the above facts in mind, it is possible to 
write the basic equations for a split ring with fixed 
ends. 


Basic Equations for a Split Ring with Fixed Ends 


The equations for X,,, X,, and X, are symbolically 

the same as Eqs. (la), (2a), and (3a), the only difference 
being in the limits of integration. The limits are now 
s and s/2, or in terms of the parameter 0, 0 = a’ to 
6= — w/2. : 
This change of limits requires that the five basic 
integrals of Step III be evaluated again. The terminol- 
ogy is the same as that of Step III with the addition of 
the following expressions: 


Ky _ f do/ V1 — k? sin? 6 
Ky = f V1 — k? sin? 6 do 


s/2 
J ds = b(Q2 + Kz) 
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(g) J yds = b2(Qs + Ke) + 
b°k LS ae 1 — k?\,, 
cos @ ¢— + cos? a+ ( a Yor x 


2 Rk? 
1 — k 
+ cos? ) — log. | | 
a y re 


Ee (cos a + in k 
k? 

















» 








i. /1 Bi 
y2 Sin ksina + Vi: —j- Re sin=* (—2&) 


s/2 
(i) I y? ds = b® (Q2. + Ke) + b*k cos a 





1 — Rk? 1 — k? 
> - + cos? a + i) 53 tos. cos a + 


1 — 1 — 
Re + cos? a] — log, a. —-{— 
3 


3 ye a — (1 + k*) sinta + Rk? sina + 


2k? — | - 1 — k? 
b' ( 3k ) (Q2 + Ke) + ( 3B ) (Q) + | 


(i) J " gtds = a°(Q2 + Ks) + 


ao. | es 
3 aie — (1+ k?) sinta — k*?sin'a — 


2k? — | 1 — kt 
ab (=) (Q2 + Ke) + ( 3 )@ + xi 


These equations are solved for various values of a, 
and the results plotted at the end of this paper. In 
applying these results, it is to be remembered that the 
origin was chosen at the top of the ellipse. Although, 
because the ring is split, the origin is now only a point 
in space, this does not preclude its being the point of 
action of the unknown loads. It is to be remembered 
that any point in space could have been chosen as 
the origin, provided that the values of the integrals 
were calculated with reference to that point. Eqs. 
(f) through (j) are derived on the basis of the origin’s 
being at the top of the hypothetic ellipse. Therefore 
the substitution of values from these equations ((f) 
through (j)) in Eqs. (la), (2a), and (3a) yields a set of 
redundants which act at the top of the hypothetic 
ellipse. This fact should be kept in mind when deter- 
mining the effects of X,,, X,, and X, on the determinate 
loads on the frame. 























STEP V 


In many cases, complicated loadings may be broken 
down into the simple loadings of Steps VI and VII, 
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and the principle of superposition applied. When this Panam i 
is not possible, the general equations of Steps II and IV +/2 
must be applied directly. Evaluation of f ds/b 
































Step VI k=0 0.1737 0. 500 0.707 0. 866 
pa a 
When the loads on both sides are equal and sym- 0 3.1416 3.1178 2.9350 2.7012 2.4222 
metrically placed, the integrals containing Mp have the 10 2.9671 2.9459 2.7836 2.5772 2.3336 
same numerical value on either side of the cut and the 20 2.7925 2.77388 = 2.6307 2.4496 = 2.2377 
left-and right-hand summations may be combined. n = oe apo + a oa “« pee 
It should be noted that under these conditions the 50 > ane one 2. 1526 20221 «1. 8686 
integrals f° Mo ds and Moy ds need only be calcu- = go. 944 = 2.0818 —:1.9854 1.8626 1.7172 
lated for one-half of the ring and that the integral 70 1.9199 1.9078 1.8148 1.6962 1.5549 
JS Me ds disappears from the equations. 80 1.7453 1.7334 =1.6418 =—-:1.5247 1.3850 
The equations of Steps II and IVb become: 90 1.5708 1.5580 = 1.4675 = 1.3806 ~— 1.111 
x = (4 fy ds)( f Moy ds) — (4 f-y? ds)(f/ Mo ds) TABLE 2 
Ag = 's/2 
D Evaluation of J y ds/b® 
(1b) St 
XY. = (4 fy ds)(f/ Mods) — (4 fds)( f’ Moy ds) k =0 0.1737 0.500 0.707 0.866 «° 
=—— a’ 
D ob 0 3.1416 3.1178 2.9350 2.7012 2.4222 90° 
(2 ) 10 3.1407 3.1069 2.9329 2.7008 2.4214 80° 
oa? (3b) 20 3.1345 3.1044 2.9323 2.6962 2.4182 70° 
— 30 3.1180 3.0832 2.9137 2.6831 2.4079 60° 
iii Da ae 2ds) — (2 f’y ds)? 40 3.0863 3.0572 2.8854 2.6575 2.3855 50° 
( S ( S9 50 3.0349 3.0046 2.8362 2.6122 2.3446 40° 
The equation of Step [Va becomes X, = 0. 60 2.9604 2.9346 2.7658 2.5444 2.2800 30° 
No attempt has been made to evaluate the integrals 7 2.8596 2.8249 2.6670 2.4489 2.1860 20° 
7 ‘fic loading condi- 80 2.7301 2.7072 2.5374 2.3209 2.0601 10° 
S'Mr ds J Mey ds for any specific loading cot 90 2.5708 2.5446 2.3796 2.1656 1.9011 0° 
tion because there are any number of load distributions 
that may be classified as ‘‘symmetric loadings,” and, TABLE 3 
as stated in the abstract, it is not the purpose of this ; 6/2 
paper to give a number of specific solutions. Evaluation of f x ds/ab 
. . ss Ss 
If the loading is such that the integrals containing My : 
do not lend themselves to simple analytical solution, k=0 0.1737 0.500 0.707 0.866 a° 
. . a 
standard numerical or graphic methods should be used. 0 2.0000 1.9897 1.9132 1.8181 1.7092 90° 
10 1.9848 1.9747 1.9003 1.8067 1.7011 80° 
Step VII 20 1.9397 1.9301 1.8608 1.7737 1.6762 70° 
, : 30 1.8660 1.8575 1.7954 1.7169 1.6303 60° 
When the loads on both sides are symmetrically 40 1.7660 1.7586 1.7037 1.6353 1.5597 50° 
placed but at in opposite directions, it is possible (as 50 1.6428 1.6362 1.5884 1.5285 1.4622 40° 
in Step VI) to combine the left- and right-hand sum- 60 1.5000 1.4940 1.4516 1.3981 1.3383 30° 
mations. In this case it should be noted that the 70 1.3420 1.3367 1.2972 1.2473 1.1913 20° 
. : 80 1.1736 1.1683 1.1301 1.0817 1.0274 10° 
integrals JM, ds and Jf May ds disappear and that => 1.0000 0.9947 0.9569 0.9088 0.8544 0° 
the integral “ Mox ds need only be calculated for one : , 
side. TABLE 4 
° 's/2 
The equations of Steps II and IVb become Desfestion of J y? ds/b3 
X, = 0 (1c) “ 
‘ k=0 0.1737 0.500 0.707 0.866 a® 
= 2 
X, = 0 (2c) * 
2 {Mx ds Me ds 0 4.7124 4.6672 4.3498 3.9375 3.4396 90° 
xX, a= af Nak — i Mew. (3c) 10 4.7124 4.6487 4.3473 3.9375 3.4389 30° 
2f'x? ds Sx? ds 20 4.7121 4.6570 4.3566 3.9874 3.4494 70° 
- 30 4.7105 4.6448 4.3475 3.9352 3.4383 60° 
The equation of Step IVa becomes 40 4.7047 4.6506 4.3433 3.9313 3.4344 50° 
2S Mex ds 50 4.6892 4.6314 4.3267 3.9162 3.4217 40° 
X, =— - 2 coceaees — = 60 4.6571 4.6082 4.2981 3.8879 3.3943 30° 
2 fx? ds + (2/3)c%d*I/L*I, 70 4.5986 4.5691 4.2404 3.8333 3.3393 20° 
S Mex ds 80 4.5021 4.4589 4.1423 3.7369 3.2458 10° 
4.3562 4.3058 3.9992 3.5986 3.0998 0° 





fx? ds + (1/3)c%1/LI, _ 
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TABLE 5 


's/2 
Evaluation of I X? ds/a%b 
St 


k=0 0.1737 0.500 0.707 0.866 a 





0 1.5708 1.5684 1.5202 1.4651 1.4048 90° 
10 1.5691 1.5651 1.5185 1.4639 1.4039 80° 
20 1.5570 1.5524 1.5080 1.4550 1.3970 70° 
30 1.5255 1.5216 1.4799 1.4310 1.3775 60° 
40 1.4680 1.4638 1.4275 1.3837 1.3366 50° 
50 1.3807 1.3778 1.3457 1.3082 1.2675 40° 
60 1.2637 1.2610 1.2335 1.2009 1.1657 30° 
70 1.1207 1.1186 1.0936 1.0645 1.0227 20° 
80 0.9582 0.9555 0.9325 0.9049 0.8744 10° 
90 0.7854 0.7834 0.7600 0.7326 0.7024 0° 





Evaluation of Integrals for Values of k and a’ 


Curves are plotted for a = 0°, 10°, 20°, 30°, 40°, 
50°, 60°, 70°, 80°, 90°. 

The value of the integral in question is taken as the 
ordinate and 6 = sin! k is taken as the abscissa. It 
should be noted at this point that 8 = O represents a 
circle. 

The integrals were evaluated for the five values of 8 
given below. 


k? k B 
0.000 0.0000 0° 
0.030 0.1737 10° 
0.250 0.5000 30° 
0.500 0.7070 45° 
0.750 0.8660 60° 


The results of the evaluations in Tables 1—5 are given 
in a dimensionless form and may be applied to an ellipse 
of any size. 


Use of Curves 

The curves given in Figs. 5, 6, 7, 8, and 9 are general 
and represent (1) the length of an elliptical arc, (2) 
the moment of length about the X axis, (3) the moment 
of length about the Y axis, (4) the moment of inertia of 
the length about the X axis, and (5) the moment of 
inertia of the length about the Y axis. 

The denominators of the expressions for X,,, X,, and 
X, are completely determined from these curves, but 
only part of the numerators can be found from 
them. 

It may be possible (depending on the type of loading) 
to evaluate the terms containing /) in a manner similar 
to that used for the terms depending only on the shape, 
but it is advisable when the loading is complicated to 
determine the values graphically. 

These curves may also be used for rings of oval shape 
when there is not too much deviation from an elliptical 
shape. 


Range of Curves 


a/b varies from 1 to '/2 and a@ varies from 0° to 
90°. 








COC 6 2 HH HW TC 60 
& sink 


Fic. 5. fds/b vs. 8 = sin! k. 





90° 50° 60 





/O 
A Sink 
Fic. 6. fy ds/b? vs. 8B = sin=! k. 


Application of Formulas to Fuselage Rings 


Before it is possible to determine the loads in the 
determinate structure, it is necessary to put the bulk- 
head in balance by skin reactions. This may easily be 
done by using the basic concepts of shear flow, keeping 
in mind that the shear distribution over a section in the 
region of a concentrated load does not exactly follow the 
law that:? 
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50’ #0" 50° 60° 


Oo 
oO fF 
B- siak 
Fic. 7. fx ds/ab vs. B = sin! k. 
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1° £0 3°” 0° KH 6 
B- sin-'k 
Fic. 9. x? ds/ab vs. 8B = sin! k. 
q = VO/I 
where: 
q = shear in pounds per inch 
V = vertical shear in pounds 
I = moment of inertia of section 
Q = static moment of area 
The shear distribution determined by the above 


formula would give a maximum running shear in the 







































































125 
4.8 la-0%4 % 203 30h 40° 
44 -— 9 Sy 50% wees 
‘re 
N 4.0 SSN 403 
8 36 . RSS 
AS 
Qs 3.2 ——1 $j} __— RS 
28 panei =e ‘E 
24 ——-+ ales 
Sl a 2 — EE 
16 —_—_}-—— a eee Oe 
b2}-— —-— +—_——+- — 
a Sees <a 
¥ —_ —— al EEE 
ra] 
Cr Tra ate aT Oe 


- sia'k 
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Fic. 10. The shear center for the assumed distribution is 
given by the formula h = a/2 + 2A/y, where h = distance from 
Y axis to shear center, A = shaded area, and y = projection of 
effective length on the Y axis. 


skin at the neutral axis of the skin cross section. Ac- 
tually, the running shear loads are nearly constant over 
the entire periphery, and it is sufficiently accurate for 
purposes of bulkhead analysis to assume constant run- 
ning load over an effective length of skin. The effec- 
tive length of skin may be assumed to be that length 
of skin outside a line drawn parallel to the Y axis 
through x = a/2 (see Fig. 10). The projected shear 
per inch is determined by dividing the load required 
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The integrals to evaluate are: fds, fy ds, fy? ds, 
JS Mo ds, and [Moy ds. 

The values of the first three integrals are found 
directly from the curves or tables. 


k? = [(2a)? — a*]/(2a)? = 0.75 
k 0.866 
8B = sin k& = 60° 


The following values of the integrals are determined 
from Tables 1, 2, and 4. 


— «/2 
ds 2.42226 
I x/2 
— «/2 
Fa y ds = 2.42225? 
+ 2/2 


— «/2 
I y? ds = 3.439603 
+ x/2 


The integrals containing W/o are evaluated as follows: 
My = Pv = Pbsin #6 
ds = bv/] — k? sin? 9 dé 
Poe [sin oT = sin 6 d0 
0 








Fic. 11. Loading for illustrative example. 


— x/2 

for balancing (applied at the shear center of the effec- f My ds 
tive skin segment) by the projection of the length of " 
the segment along the Y axis. The actual shear per Performing the above integration, 
inch in the skin at any point is equal to the projected 
shear per inch as defined above. 

Once the ring has been balanced, the determinate 
moments, shears, and axial loads (Wo, Vo, and Po) may 


—x/2 
f My ds = Pb*(0.69006) 
0 


y = d(1 — sin 8) 





be determined at any point by treating each half of the —x/2 —1/2 
ring as a cantilever with its free end at the cut. J Moy ds = Pos sin 0(1 — sin @) X 
All the integrals may now be evaluated and X,,, X,, . : a 
' : V1 — k sin? 6 de 


and X, determined. The total loads at any point are 


found to be the algebraic addition of the effects of fo a/2 ——___ 
Xm Xp and X, to Mo, Vo, and Py. = Pb* f sin 64/1 — k* sin? 6d6 — 


—x/2 Le 
ILLUSTRATIVE EXAMPLE J sin? 6+/ 1 — k? sin? 6 dé 
0 


The problem shown in Fig. 11 is now solved to illus- _,, , , , 
‘ = Performing the above integration, 
trate the use of the tables and curves. Since the ring : 
and loading are symmetric, the equations developed f a fy ds = 1.19880PB 
under Step VI are used. 0 








Substituting the above results in Eqs. (1b), (2b), and (3b), 








Xp — 412-4228") (11988075) — 4(3.439604)(0.69006P5*) _ 9 915155 
4(2.4222b) (3.4396b°) — [2(2.4222b?)]? 

a 4(2.4222b?)(0.69006Pb2) — 4(2.4222b) (1.19880) Pb? = 0.50004P 
4( 2.42225) (3.439608) — [2(2.4222b%)]? 

Xp = 0 


It should be noted that Xp could have been determined by inspection as equal to P/2, due to the fact that the 
loading is symmetric about two axes. 
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Vision, Hearing, and Aeronautical Design 
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C IS PRINCIPALLY through the sense of vision that 
pilots are aware of the nature of terrain, weather 
conditions, the terrestrial orientation of their aircraft, 
and the position and movements of enemy planes. 
Even so-called ‘‘blind flying”’ is primarily a visual func- 
tion. Here the “contact” is transferred from the 
horizon and ground objects to the special instruments 
upon the panel. Hearing likewise is of critical im- 
portance to military flyers in that command informa- 
tion and assistance from ground control stations must 
be audible over the communication system. Also, 
intercommunication systems within aircraft must 
function as the principal means of maintaining prompt 
and coordinated action of members of air crews in 
carrying out military missions effectively. These 
facts are generally known and accepted by plane de- 
signers. However, engineering success in providing 
for ease of seeing and ease of hearing is quite variable, 
and these basic pilot requirements merit continued 
attention. 

Take, for example, the almost complete blanking off 
from vision of an area directly ahead of the pilot when 
a plane changes its normal flight attitude to a landing 
attitude. This is particularly vicious in the case of 
fighter aircraft with large air-cooled engines in the nose 
of the ship, and the condition has resulted in a great 
number of avoidable landing, taxiing, and take-off 
accidents. Aircraft carrier accidents have undoubt- 
edly resulted from this blanked-off area, in spite of the 
fact that a landing signal officer is supposed to direct 
the landing properly. Mirror reflections of the area 
immediately ahead of the plane could be provided. 
Tricycle landing gear on some few aircraft have done 
much to benefit this blindness due to radical change in 
plane attitude. 

Engineers are no doubt aware of the desirability of 
placing pilots, observers, and gunners as near as possible 
to the transparent panels through which they must get 
their visual impressions. The nearer the eye is to the 
window, the larger the visual field that can be seen 
through it or, for a given angular field, the smaller the 
window needed. The ordinary spectacle lens 1'/2 in. 
in diameter gives the wearer an uninterrupted view of 
almost 90°. Proximity to airplane windows also 
makes for clear visual fields: (1) It is easier to provide 
the requisite transparent area free from obstructing 
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(2) 
marks, specks, and scratches on the window surfaces 
are less bothersome because they are out of focus; 
and (3) when near the window, the observer’s head 
makes an area of shadow which is helpful, particularly 
at night, since it blocks off reflections from the inner 
surface. 


pillars, panel joints, and opaque accessories; 


THE GUNNER'S VISUAL PROBLEM 


The position of gunners in relation to windows is 
especially worthy of consideration. Take a hypo- 
thetic set of measurements for discussion. The 
gunner’s eye must be about 9 in. behind the gunsight; 
the sight is 20 in. from the aiming panel; the gunner’s 
seat is therefore placed against the back wall of the 
turret, and the space in front of him is well filled with 
gun supporting and operating structures. The aiming 
panel is a beautifully clear, optically perfect surface 
14 by 14in. The engineers are to be congratulated on 
this window but not on the position of the gunner. 
For him the situation resembles that met in tunnel 
vision. His eyes are 29 in. away from this superb 
window which can provide him with a view of only 
28°—that is, 28° of the full 60° to 90° that are so 
critically needed. 

The aiming panel must, of course, be well mounted, 
and this calls for a rather heavy metal frame. Some- 
times it is as much as 2 in. wide. Other small, usually 
curved windows border the aiming panel, and each of 
these must have or share a frame with adjoining win- 
dows or wall structures. When the guns are in front 
of the gunner, sections of the mounting structure, 
brackets, electrical switches, etc., frequently have 
locations that obscure portions of the windows from 
the gunner’s eyes. The total effect is a lattice window 
with the lattice predominating over the clear space in 
some instances. From the gunner’s view through a 
particular panel, a total of only 25 per cent is un- 
obstructed to binocular vision; 53 per cent is totally 
obstructed; 22 per cent is obstructed for one eye. 
For the gunner the irregular lattice is a hazard to vision 
and life. It is true enough that through the lattice he 
can see a lot of landscape as he looks from one opening 
to another. What he needs to see, however, is moving 
planes. When his eyes catch sight of one, they try to 
follow it. This is done by a special kind of coordi- 
nation in the action of eye muscles and is known as 
pursuit eye movement. The eye actively glides along 
at a rate set by the visual target and appropriate to 
keep it in clearest possible view on the fovea. Visibil- 
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ity of the moving target is essential for this type of 
seeing. Whenever the moving target momentarily 
disappears, as happens when it passes behind opaque 
lattice areas, the pursuit movement is immediately 
stopped. The eye jumps to the far edge of the ob- 
struction and waits. When the target emerges, a new 
pursuit movement must be organized by the nervous 
system. This necessarily costs an average delay of 
0.15 of a precious second before clear pursuit vision 
can be re-established. The gunner is, of course, not 
completely blind during this interval, but his visual 
efficiency on the moving target is hampered by every 
opaque gap across which he is compelled to operate. 

In turrets where the guns may be mounted at the 
sides and rather low down, it should be possible to 
reduce the amount of structure in front of the gunner 
and to bring him forward much nearer the front panel. 
This will increase his angle of uninterrupted view, 
make visual pursuit of his targets easier, and reduce 
the blinding effect from the flashes of his own guns. 


VISION INSIDE THE COCKPIT 


Inside the cockpit there is a confusing multiplicity 
of instruments, and a lack of a suitable grouping of 
essential flying instruments according to a well-planned 
pattern for visual perception is noticed. For purposes 
of night flying, where visual function is largely de- 
pendent upon maintenance of dark adaptation of the 
eye, instrument panels in general are subject to the 
following serious faults: 

1. There is too large an illuminated area. Instru- 
ments that are referred to infrequently are constantly 
illuminated. 

2. The color of the transmitter or reflected light 
from these instruments is usually not of the spectral 
band least disturbing to night vision. 

3. Intensity, as well as total area of illumination, 
is considerably too high. This is usually a criticism 
equally applicable to direct or indirect lighting, radio- 
luminous, or fluorescent marked instruments. The 
only light that can be controlled easily both as to 
spectral character and brightness level is indirect red 
light. The use of a spectral red whose transmission 
band lies in the region of 600 millimicrons is recom- 
mended for instrument lighting... The intensity of such 
light may be permitted to vary considerably without 
adverse effect upon night vision. 

Transparent cockpit enclosures, although of fairly 
clear optical plastic, discolor with explosure to sun- 
light; some of them become finely checked because of 
temperature changes and vibration. Nearly all of 
them are capable of being scratched too easily, and the 
result is interference with vision. Measurements of 
visibility through plastic windows, compared with 
bulletproof glass and with open cockpit yiew, show that 
the glass produces a slight loss and that the loss from 
the plastic may be 5 or 6 times as much. Flat panes 
consistently provide better visibility than do curved 
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plastic surfaces, but they may also give an increased 
drag. A compromise must be made between visibility 
and aerodynamic considerations. Rapid strides in 
development of better plastics are being made, and it 
should be possible soon to mold transparent cockpit 
enclosures of better grades of optical plastics in one 
piece. Surface hardening of such molded parts is 
desirable. Close attention should be given to reducing 
inside reflections that are particularly troublesome 
from concave surfaces. 

In military aircraft it is often found that vision for- 
ward and to either side is fairly well provided for but 
that, even though the pilot may have the usual well- 
developed ‘“‘rubber neck,”’ he is handicapped by certain 
design features of his plane in seeing what may be above 
or aft. It is often forgotten that in actual combat 
pilots refuse to leave the “‘greenhouse’”’ closed and must 
therefore use goggles. Any fixed aperture made 
available for the pilot’s face should not be too small to 
permit use by a bespectacled or goggled aviator. 


AIRSICKNESS AND VISION 


Airsickness of passengers and crew may result from 
unavoidable motion stimulation of the vestibular 
mechanism of the ear; from rapidly changing gravita- 
tional forces acting on viscera, muscles, and joints; and 
from apprehension and past unhappy memories of plane 
travel. All of these upsetting stimuli are, as a rule, less 
disturbing if those affected can see out and establish vis- 
ual contact with the horizon, with cloud formations, and 
with the ground scene below. The old instructions to 
novices in flying should be remembered: ‘‘Never look 
at the up wing; watch the down wing.’’ Many troop- 
carrying glider craft afford virtually no opportunity 
to look out and establish visual contact beyond the 
plane. To arrive at the scene of battle with a load of 
thoroughly ill troops contributes nothing to fighting 
morale and effectiveness. A part of the answer to this 
difficulty is—don’t require passengers and crew to fly 
blind. 


HEARING SUFFERS MoreE THAN SIGHT 


The ears more than the eyes are subjected to en- 
vironmental stress through flying. It seems to be 
true of modern aviation that every time the engineers 
increase the power and speed of our airplanes, the ears 
of the pilots take a greater beating. Although the ear 
is a magnificent little mechanism—the most intricate 
mechanical structure in the human body—1tt is a rather 
delicate device and one that seems ill designed for 
modern war. But the ear has gone to war along with 
the rest of the soldier, and the service it renders in the 
face of acoustic stress is to be admired. 

Airplanes have always been noisy, and they are 
becoming noisier. A thousand horsepower fed into a 
propeller is able to agitate the atmosphere in a thunder- 
ous manner, and when the engine delivers 2,000 hp. the 
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din is doubled, or actually more than doubled, because 
as the tip speed of the propeller increases a larger pro- 
portion of its driving energy is converted into sound. 
When this energy pounds on the ear, it is striking a 
mechanism so sensitive that less than one quintillionth 
of a horsepower is needed to produce a faint sensation 
of hearing. In addition, more horsepower means more 
speed and hence, more turbulence about the ship. It 
is this turbulence of the slipstream over the wings and 
about the fuselage which produces the distressing, high- 
frequency random noises that sound like a mighty 
“shhhhhh.” In some respects the noise from the 
turbulence about the plane is more of a problem than 
is the noise from the propeller itself. This is demon- 
strated in planes that do not have propellers. Con- 
trary to popular notions, the interior of a glider plane 
moving at about 150 miles per hour is a noisy place. 
The noise level is about 115 db., and conversation in 
such a place is difficult if not impossible. 

In any really fast-moving vehicle the noise is ran- 
dom; that is to say, all frequencies of the spectrum are 
present. When one listens on the ground to a plane 
high overhead one hears only the low frequencies of the 
propeller. But inside the plane it is different; there 
all frequencies added together at once are heard, pro- 
ducing a noise that is to sound what white light is to 
light. And as a general rule, the greater the speed, the 
“whiter” the noise. Also, as a general rule, the whiter 
the noisc, the more objectionable it is to the ear. 
“White” noise is objectionable for three reasons: 
(1) It is disagreeable; (2) it produces temporary deaf- 
ness; and (3) it spoils communications. 

That ‘‘white” noise is annoying needs little argu- 
ment. No one has been found who really enjoys it. 
It is true, however, that attempts to prove that long 
exposure to intense airplane noise is damaging to 
human efficiency have produced essentially negative 
results. When adequately motivated, a man can code 
a message, add columns of figures, coordinate his move- 
ments, react to a signal, etc., about as well after 8 
hours in a noise of 115 db. as after a similar period in 
the quiet. Despite this remarkable experimental re- 
sult, the subjects all report that they find the noise 
unpleasant, and they are happy when it is turned off. 


TEMPORARY DEAFNESS 


That airplane noise produces deafness is a well- 
known fact. In normal ears this deafness shows two 
characteristics: It is restricted more or less to the high 
frequencies and, fortunately, it is usually temporary. 
After 8 hours in an airplane noise of 115 db., the normal 
ear shows a hearing loss of about 40 db. in the region of 
4,000 cycles. It has sometimes been supposed that 
low-frequency airplane noise produces high-frequency 
hearing loss. On the contrary, it is the high-frequency 
components of the noise that produce the high-fre- 
quency loss. The ear, for some unknown reason, is 
more vulnerable at these high frequencies. 
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Recovery from a 40-db. hearing loss usually occurs 
in about 24 hours. The recovery is rapid at first and 
then proceeds more slowly. About half of the loss is 
regained by a normal ear in the first 3 hours after ex- 
posure. Some ears apparently do not have this 
power of recovery, and repeated exposure to noise 
leaves them permanently deafened. Although - the 
noise is usually blamed for such deafness, there is a 
question as to whether the ear itself is not the real cul- 
prit. At any rate, many ears appear to develop high- 
frequency deafness without excessive exposure to 
intense noise. 


IMPROVEMENT OF COMMUNICATIONS 


The most important practical effect of airplane noise 
is the masking of communications. Not only is con- 
versation impossible in some planes, but even over 
radio and interphone speech signals are often masked 
beyond recognition. Articulation tests have shown 
that with much of the standard military interphone 
equipment a listener is able to understand only about 
50 per cent of the words spoken in the presence of an 
airplane noise of 120 db. Over the same interphones 
more than 95 per cent of the words are understood 
when there is no ambient noise present to interfere with 
the speech. 


The difficulty of communicating under the handicap 
of airplane noise calls for vigorous remedies. Con- 
structive measures can be applied in three general 
directions: 


First, the plane can be quieted to some extent either 
by improved aeronautical design or by the application 
of sound absorbent materials. Acoustic treatment 
that is light enough to be tolerated in a plane does not 
appreciably reduce the overall noise level. It does, 
however, change the spectrum of the noise by reducing 
the intensity of the high-frequency components. 
Hence, the noise in an acoustically treated plane is 
less ‘‘white,’”’ and therefore less bothersome, than the 
noise in an untreated plane. Tests have shown that 
for the same overall sound intensity conversation 
person-to-person is relatively easy in a treated plane 
but quite impossible in an untreated one. 


The second remedy calls for an improvement in the 
response characteristics of the communication equip- 
ment itself, especially of the microphones and ear- 
phones. A loud noise does not interfere with intelligi- 
bility nearly so much when instruments of high fidelity 
are used. But with microphones that have nonlinear 
distortion and with earphones that at some frequencies 
are sharply resonant, the effect of an airplane noise of 
120 db. is to reduce the intelligibility of speech by 30 
to 40 per cent. High-fidelity equipment is not yet 
being widely used in airplanes, although a few major 
improvements are now in process. Complete overall 
high fidelity from microphone to earphone must be 
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achieved if speech is to be transmitted to and from our 
most modern airplanes. 

The third remedy called for by the noise problem is 
the shielding of the microphone and the earphones 
from the noise. The oxygen mask could be so designed 
as to shield the microphone from the ambient racket, 
but many otherwise excellent masks suffer from acous- 
tic defects. This problem is now under study, and 
improved noise shields for hand-held microphones are 
being developed. The earphones and the ear of the 
listener can be shielded from the airplane noise by 
means of an acoustic socket designed to provide a 
tight seal against the side of the head. In present 
military equipment this provision has been neglected, 
but improved devices are now in production. Some of 
these ear protectors serve to reduce the unwanted sound 
in the aviator’s ear by 20 to 50 db. 


CONCLUSIONS 


In general, it can be said that the problems raised by 
intense ambient noise are serious but not insoluble. 
Judicious use of sound treatment in the plane, con- 
version to high-fidelity microphones and earphones, 
and the development of acoustic devices to shield the 
mouth and the ears-of the personnel will permit the 
aviator to carry on in the best noises that the aero- 
nautical engineers are now planning to produce. 


(Continued from page 114) 
F = [1/(2 sin y)][sin (2 + l)y + sin ny — sin y] 


With the aid of this formula the two sums can be calculated. 
Substituting m for n and 27/(m + 1) for y, 
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It can be shown in a similar way that 


m 


cos [4ri/(m + 1)] = —1 


i=1 
Consequently, 
N = (8/8)(m + 1) 
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The flying and efficient fighting of modern planes are 
largely dependent on the special sense of vision. The 
eyes fortunately suffer no great decrement in function 
from the swift movement or high and changing ele- 
vation of the airplane. Vision is adequate to the 
basic task assigned to flying personnel. The chief 
difficulty is in providing for optimal visibility through 
the structures of the plane and for continuous visual 
check on the environment surrounding it in both day 
and night flying. Ideal visual conditions are not wholly 
attainable because of aerodynamic and structural neces- 
sities. However, the effort should be made to gain all 
possible visual advantages. The problem is a continu- 
ing one and offers important strategic possibilities. 

Both seeing and hearing, if accompanied by pro- 
longed attentive effort, especially under conditions of 
unfavorable plane design, are capable of contributing 
fatigue and loss of efficiency to pilot and air crew. It 
has been proved worth while to give the airplane engine 
an adequate combustive mixture by supercharging and 
to pay special attention to protecting the oil in the en- 
gine against “‘foaming”’ at high altitudes and reduced 
barometric pressures. It is proving and will prove 
worth while to consider the flier’s eyes and ears and the 
rest of his mortal body and to reduce in every possible 
way the tendency to physiologic and psychologic 
“foaming” in him. 
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A Method for Rapid Estimation of 


Helicopter Performance 
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SUMMARY 


Several papers are available which treat the analysis of rotating 
wing aircraft with considerable rigor and completeness. While 
these treatments are satisfactory for exhaustive analyses of 
specific helicopters, they are exceedingly complex. Since simple 
methods of performance estimation are available for airplanes, 
it is believed that a comparable system for helicopters should 
prove useful for comparative performance studies, preliminary 
estimation of characteristics of proposed designs, and determina- 
tion of the effects of small changes in design. Frequently, only 
limited data are available for studies of this type, making elabo- 
rate analyses impossible or undesirable. The present report pro- 
poses a system of performance estimation which, it is hoped, will 
fulfill the need. The method is relatively simple and rapid. 

The method is based on the conception of the rotor as an ac- 
tuator disc and the well-known momentum relations for an ac- 
tuator disc. All of the major variables in design and operation 
are taken into account. Comparisons have indicated that the 
results obtained by this system are good approximations to the 
results of more elaborate analysis. 

In many respects this paper was suggested by Knight’s paper,’ 
but the method of approach is different and somewhat greater 
generality is obtained. 


SyMBOLsS USED 


B- = tip loss factor, assumed equal to 0.95 

6 = number of blades 

Cr = thrust coefficient, Cr = T/(pw?rR*/2) 

c¢ = rotor blade chord at r/R = 0.7 

f = equivalent parasite drag area, sq_ft. 

K = ratio of effective power available at main rotor shaft 
to total brake power 

P = engine brake power, ft.lbs. per sec. 

P; = power available to produce thrust 

Pr = power required for rotor torque due to profile drag 

R= rotor radius, ft. 

S = resultant axial velocity ratio (V sin a’ + u)/uo 

T = rotor normal thrust, Ibs. 

u = induced velocity at rotor, ft. per sec. 

uo = induced velocity for the hovering condition 

V = velocity of machine along its flight path, ft. per sec. 

W = gross weight, lbs. 


a = flight path angle, degrees 

a’ = angle between rotor tip plane and flight path 

6 = rotor blade profile drag coefficient at the mean effective 
lift coefficient 


= advance ratio, = V cos a’/wR 


7 

w = rotor angular velocity, rad. per sec. 

« = rotor solidity ratio, ¢ = bc/rR 

y = angle between the rotor tip plane and a horizontal 
plane 

y = blade azimuth angle (measured from down-wind posi- 


tion in direction of rotation) 








Presented at the Rotating Wing Aircraft Session, Eleventh 
Annual Meeting, I.Ae.S., New York, January 25-29, 1943. 
* Aeronautical Engineer, Research Division. 


MOMENTUM RELATIONS 


q 
| ‘HE GENERAL MOMENTUM EQUATION for the thrust of 
an “actuator disc’’» ° is 


T = 2pAuV’ (1) 


where A is the effective disc area; u, the induced "ve- 
locity; and V’, the resultant relative wind velocity at 
the disc. Referring to the helicopter rotor as shown 
in Fig. 1 

V? = (u + Vsin a’)? + V? cos? a’ 


and A = 7B?*R?,. When V = 0 
uo? = T/2prB?R? (2) 
+ 
R 
ww HORIZON _ 





UL 





Ww 


Fic. 1. Rotor geometry. 


Substituting the value of 7 from Eq. (2) and the value 
of V’ in Eq. (1) 
ut + 2u*V sin a’ + u?V? = 14 
eliminating sin a’ by the use of the following definition 
S = (V/u) sin a’ + (u/uo) (3) 


and solving for V/u, the perfectly general momentum 
relation is obtained as 


V/uo = VW (uo/u)? + (u/uo)? — 2(u/U9)S (4) 


This form of the momentum equation is plotted in Fig. 
2, which is used in estimating performance. For reas- 
onable values of S, when u/u < 0.20, Eq. (4) reduces to 





V/Uo = to/u 


In hovering flight, 7 = W. Then from Eq. (2) 
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Fic. 2. Solution of the momentum equation for helicopters (Eq. 4). 
V/uo 
W \' t/a 
u = 15.26 (=r) (o»/») (5) 


It will become evident that S is a useful parameter. 


RoToOR SPEED 


Since y is always a small angle, 7 = W is a good 
approximation for all flight conditions. Then 





ae a 
o prw*R4a/2 
and 
2 Va 
o= ts: A (6) 
prR‘a(Cr/o) 
This can be written 
16.36 (W Pa y” 1 
pie (“) pole) (Cr 7a) 7) 


Sissingh’ has shown that the static efficiency of a 
rotor is a maximum when C;/o = 0.20. This relation is 
very nearly true for any reasonable rotor design and is 
roughly true even at moderate forward speeds. Marked 
improvement in airfoil characteristics may result in 
some change in the optimum thrust coefficient, but the 
relation is sufficiently accurate to determine the opti- 
mum rotor speed at low rates of advance at any given 
flight condition. The rotor speed will be determined by 
this relation at some flight condition, probably hovering 
at sea level, and the rotor r.p.m. at other flight condi- 
tions will probably follow from considerations ot engine 
operation. Substituting the value of C;/c, the opti- 


— 


Po 


_ 36.6 


= (8) 


w 


p 


VAL, 


V/uo = V(uo/u)?+(u/uo)? —2(u/uo)S, when u/uo < 0.20 


= Uo/u. 


PROFILE TORQUE 


The profile torque of a rotor is 


2r R 
Or = 2. fi J 5(p/2)ber dr (wr + U sin y)*dy 
2r Jo 0 


where U is the velocity in the plane of the rotor. U 
Vcos a’. Integrating 


Or = (p/2)(8/4) omw*R® (1 + pu?) 
The power absorbed by the profile drag is 
Pp = Qe = (p/2)(6/4) omw®R> (1 + u*) = (9) 


Sissingh® 7 has shown that the mean lift coefficient is 
approximately 


C, = (3/B*)(Cr/o) 


from which 


6 in Eq. (9) is the profile drag coefficient of a rotor blade 
section at the section lift coefficient obtained from Eq. 
(10). This relation for the mean lift coefficient is, to be 
sure, only an approximation and is derived for the case 
of static lift, but it provides a means of estimating the 
effective blade profile drag coefficient. For this purpose 
it is sufficiently accurate and tends to be conservative 
for flight at moderate speeds. 


POWER AVAILABLE FOR THRUST 


The power left for producing thrust after the losses 
are accounted for is 


P, (11) 


550K b.hp. — Pr 











ESTIMATION OF HELICOPTER PERFORMANCE 


The power P; is the power available to the rotor to 
produce a jet of air. The power imparted to the jet is 
equal to the thrust times the total velocity in the di- 
rection of the thrust. 


P, = T(u + V sin a’) 
and since 7 = W 
P, = Wu + WV sin a’ (12) 
Equating Eqs. (11) and (12) and dividing by Wu, 
Uu V 550K b.hp. Pr 
—+ -si = —_————_ - — 13 
Ug Ug —_ Wuo Wuo (13) 


The left-hand member of this equation is, from Eq. 
(3), equal to S. Then 
_ 550K bhp. Pr 
Wuo Wu 
Substituting the value of Pz from Eq. (9), putting the 


value of w from Eq. (6) in the expression for Pz, and sub- 
stituting the value of m from Eq. (2), 


550K bhp. 8B 
Wug 2(Cr/o)¥? 


S (14) 


S= 





7 (1+ 4%) (15) 


Eq. (15) is the useful equation for performance esti- 
mation. From these equations it is evident that the 
useful parameter S has two meanings. It is the ratio 
of the total axial velocity component to the axial 
velocity when V = 0, and it is also a measure of the 
power available for the production of thrust after all 
losses are accounted for. It appears in the latter sense 
in Eq. (15), which can be easily evaluated. The quan- 
tity 4 may not be known at the beginning of a specific 
performance calculation, but its effect is small and an 
assumed value is quite satisfactory. Having calcu- 
lated a value of S, Eq. (4) or Fig. 2 representing Eq. 
(4) gives a series of corresponding values of V/um and 
u/u. These values define a series of flight conditions. 
Then the angle a’ is defined by the equation 


S = (V/u) sin a’ + (u/s) (16) 


FLIGHT PATH 


The flight path is, from Fig. 1, 


ER RIE Se Ts, TS OR (17) 
i WwW 
where D is the total drag of the machine. D = Dz + 


Dp. 
The horizontal component of the parasite drag of the 
machine is 


D, = f(p/2)U* = f(p/2)(V cos a’)? 
The drag due to the profile drag of the rotor blades is 


2 R 
Dy a 3 ¥ I 8(p/2)bcdr (Usiny + wr)? sin y dy 
2r J0 0 
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(p/2)(6/2)bc UwR? 
= (p/2)(5/2)orUwR?® 
Adding Dz and Dp and dividing by W 


_ pu [6 


_ pV cosa’ (6 , ee , 
= PV 008 0 (3 ona + fl cos a’) (18) 


The flight path angle is then, from Eqs. (17) and 
(18), 


a= 


, _ pV cosa’ (5 . s : 


The angles are in radians. The equation appearing in 
the summary gives the angle in degrees. 
VERTICAL CLIMB 
In vertical climb, Eq. (3) reduces to 
V u u V 
=§ 2 e's @ 2 wS-< (20) 
Uo Uo Uo io 


Writing Eq. (4) in the form 


(J (c+ (-2()s 
Uo u Uo Uo 


the value of u/u can be substituted from Eq. (20). 
Solving for V/u in terms of S, the solution is V/u = 
S — (1/S) or V = um |S — (1/S)]. 


Since V is a vertical velocity, the vertical rate of 
climb in feet per minute is 


Veg = 60u[S — (1/S)] 
Ce » 


(21) 


the value being taken for the condition 
w= 0. 
From this equation, it is apparent that the hovering 


condition is S = 1.00. 


EVALUATION OF K 


The power losses in per cent of total b.hp. are ap- 
proximately as follows: 


Transmission losses 31/2 per cent 
Cooling (blower) power 4'/, (unsupercharged) 
Slipstream rotation 5 


Tail rotor 6 


Total loss 19 per cent 


For supercharged engines, cooling power may be about 
3 per cent at sea level. 


K = 1.00 — 0.19 = 0.81 


The most convenient procedure is to estimate the 
value of K for one flight condition and use the same 
value for all conditions. A method of estimating the 
tail rotor power loss follows. 
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ESTIMATION OF TAIL ROTOR POWER 


(For helicopters with torque counteracting auxiliary 
rotors.) The main rotor torque is 
Q = P/w = 550K b.hp./w (22) 
The thrust required at the control rotor to balance 
the main rotor torque is 


T; = Q/l 


where / is the moment arm of the tail rotor about the 
main rotor axis. The power required by an airscrew to 
deliver a given static thrust is 


550ns b.hp. = T r 
2p7R? 


(23) 


(24) 





where 7s is the static thrust efficiency, which may be 
taken as equal to 0.70. Substituting Eqs. (22) and (23) 
in Eq. (24), 


5 3/2 
Mine (= bat) 1 











wl 550 X 0.7 V2prRz? 
Simplifying, 
548 (K b.hp.\*/ aie - 
b.hp., = — ( e:) (/o) (25) 
Dr wl 


where D7, is the diameter of the control rotor and w is 
the speed of the main rotor. The diameter and rota- 
tional speed of the tail rotor should be adjusted so that 
C,/o is approximately equal to 0.20 for the hovering 
condition. K in Eq. (25) is, of course, dependent on 
b.hp.7, but an estimated K satisfies Eq. (25) with suf- 
ficient accuracy to obtain b.hp.7 and the resulting new 
and more accurate value of K. The loss of power due 
to-the tail rotor in per cent is b.hp.7/b.hp. Then 


K = 1.00 — ®& losses 


PROCEDURE 


The equations as presented do not give maximum 
speeds and rate of climb directly but give the speed 
along the flight path and the angle of the flight path 
for a series of flight stages. The determination of 
maximum rate of climb, maximum horizontal speed, 
etc., necessitates plotting a speed polar (V, vs. V,) 
from which all phases of performance at a given alti- 
tude and power can be read directly. 

The procedure in calculating a point on the perform- 
ance curve is as follows: 


1. Solve Eq. (5) for uo. 

2. Knowing the rotor speed w or using the approxi- 
mately optimum value of Eq. (8), calculate the blade 
loading coefficient C;/c. 

3. From the polar of the rotor airfoil section, ob- 
tain the blade profile drag coefficient 6 = Cp, at the 
average blade lift coefficient (Eq. (10)). 

4. Compute S from Eq. (15). 
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The value of uw is as yet unknown, but its effect is 
small and with a little practice its value can be accur- 
ately estimated for each point on the performance 
curve. When the rotor speed is that prescribed by Eq. 
(8), Eq. (15) reduces to 


_ 550K b.hp. _ 5.31 . 
Wu Vo 
5. Referring to Fig. 2, choose corresponding values 


of u/up and V/u at the computed value of S. Deter- 


mine V. 
6. Check the assumed value of yu. 


S (1 +h ?) 





uw = Vos a’/wR 


7. Then use the value of sin a’ from Eq. (16). 
8. The flight path angle is (see Eq. (19)) 


, V cos = ( 





(p0/p)W 


9. Rate of climb in feet per minute is 


V, = 60V sin a 


; orwR*® + fV cos a’) 


a=a’'—0. 


10. The corresponding horizontal component of 
speed in m.p.h. is 


V, = 0.682V cos a 


EXAMPLE 


To illustrate the results of the performance equa- 
tions, the performance of a representative helicopter 
has been calculated and plotted in Fig. 3. The climb 
curve for the corresponding airplane with both a fixed 
and a variable pitch propeller has also been calcu- 
lated. The corresponding airplane has been taken as 
one having the same brake horsepower, same span, same 
weight, and comparable care in minimizing parasite 
drag. (The parasite drag of the helicopter is taken 
equal to 1.25 times the difference between the total 
parasite drag of the airplane and the airplane wing 
profile drag. Cp, = 0.010.) 












8 


8 





V. , FT/MIN. 





V, MPH 


Fic. 3. Typical helicopter performance curve and performance 
of the corresponding airplane, both at sea level. 


Curves have been calculated both with and without 
the wu? term included in the value of S in order to illus- 
trate the magnitude of its effect. It is evident that 
for approximate work it may be neglected. 

















ESTIMATION OF HELI 


Characteristics of the machines used as examples: 


Airplane Helicopter 
WwW = 1,500 Ibs. Ww = 1,500 lbs. 
b = 33.6 ft. D = 33.6 ft. 
b.hp. = 100 b.hp. = 100 
W/S = 10.0 lbs. per sq.ft. 
Cp = 0.030 
f = 4.50 sq.ft. 
(Residual drag fr = 3.00sq.ft.) f = 1.25 X 3.00 = 3.75 sq/ft. 
‘mes. = 0.83 « = 0.070 K = 0.81 
6 = 0.0125 Cr/o = 0.20 
REFERENCES 


1 Glauert, H., Horizontal Flight of a Helicopter, R. & M. No. 
1157. 

2? Hohenemser, K., Performance of Rotating Wing Aircraft, 
N.A.C.A. T.M. No. 871, 1938. 


COPTER PERFORMANCE 135 


3 Knight, M., General Flight Analysis of the Helicopter, pre- 
sented at the Rotating Wing Aircraft Session, Tenth Annual 
Meeting, I.Ae.S., New York, January 29, 1942. 

4 Knight, M., and Hefner, R. A., Static Thrust Analysis of the 
Lifting Airscrew, N.A.C.A. T.N. No. 626, 1937. 

5 Kiissner, H. G., Helicopter Problems, N.A.C.A. T.M. No. 
827, 1937. : 

5 Sissingh, G., Contribution to the Aerodynamics of Rotating 
Wing Aircraft, N.A.C.A. T.M. No. 921, 1939. 

7 Sissingh, G., Contribution to the Aerodynamics of Rotating 
Wing Aircraft, Part IT, N.A.C.A. T.M. No. 990, 1941. 


8 Squire, H. B., The Flight of a Helicopter, R. & M. No. 1730, 
1936. 

® Wheatley, J. B., An Aerodynamic Analysis of the Autogiro 
Rotor with Comparison Between Calculated and Experimental 
Results, N.A.C.A. T R. No. 487, 1934. 


Letter to the Editor 


Dear Sir: 

In a letter to the editor from Mr. Manley on “Pendulum- 
Type Vibration Absorbers” (Journal of the Aeronautical Sciences, 
January, 1943), he pointed out the incorrectness of the assumption 
of the “‘reduced system”’ made in the writer’s paper in the Jour- 
NAL, July, 1942. He corrected this by adding one more equation 
of motion for the propeller and showed the deviation of my re- 
sult from his correct one. 

However, two further points have to be considered: One is the 
damping effect of air on propeller, and the other is the assumption, 
mk? << Iz. The latter assumption has been used in the writer’s 
paper and also employed by Mr. Manley in his work. This as- 


sumption is evidently not true, because my first approximation 
already indicates a value of » lying between 0.4 and 0.5 and his 
considerations also indicate that 4 should be made as large as 
possible. From this point of view, Mr. Manley’s calculation is 
also unsatisfactory. ‘ 

Regarding the damping effect of air on propeller, it could be 
taken into account by adding a term —c6, on the right side of the 
third equation of (1) in Mr. Manley’s letter, where c is the damp- 
ing constant depending upon the propeller speed (see Den Har- 
tog’s Mechanical Vibrations, 2nd Edition, p. 251; McGraw-Hill 
Book Co., Inc., New York, 1940). Then the amplitude of J; can 
be expressed as 





K*(K — hp?) 











Mo _ Ss nae _ mL(R + L)*p4 _ 2 Kep 2 
thas yl" al ictal Ro? — Lp? (K — Lp)? + cp? | + | (KK — 1p)? + cp 


By putting it into nondimensional form, the following equation is obtained: 





n? 1 — g? 





Ost. 
“a = yl: — n(1 + r2)g? — ni(1 — 12)8? 2 


where 
oe R a gas mk® a I 
a ete ies 
A gt 
pay Q% a ar” 
K K Mo 
2 =— 2%, = — Ou. = — 
» c : qh 7 K 


“|2 f 7]2 
— nt (1 yt +la Re 


The above equations may give a better understanding of the en- 
gine-propeller vibration characteristics. 

In addition, there is one further remark to be made. The above 
equations are valid only for the case of direct-drive propellers. 
Modern high-performance airplanes are practically all equipped 
with geared propellers. Theoretic investigation of this case in- 
volves no new principles, except more equations have to be 
handled. The calculations will be quite lengthy and cumber- 
some. 

A. H. SHIEH 
Pasadena, Calif. 
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